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INTRODUCTION 


This  is  the  final  report  on  the  research  program  on  "Radar  Target 
Di scriminat ion  and  Identi  ficat ion  using  Extinction-Pulses  and  Single-Mode 
Extraction  Signals"  supported  by  the  Naval  Air  Systems  Command  under  Contract 
N00019-85-C-0411,  covering  the  period  of  September  11,  1985  to  March  30,  1987. 
The  original  contract  was  for  the  period  of  September  11,  1985  to  September 
10,  1986,  but  it  was  extended  on  no-cost  basis  from  September  11,  1986  to 
March  30,  1987  in  order  to  complete  some  additional  new  research  topics. 

The  purpose  of  this  research  ; s  to  develop  a  new  scheme  of  ,'udar  target 
discrimination  and  identification.  This  scheme  consists  of  synthesizing 
aspect- independent  discriminant  signals,  called  extinction-pulses  (E-pulses) 
and  single-mode  extraction  signals,  which  can  be  used  to  produce  zero  response 
or  single -mode  responses  of  the  target  in  the  late -time  period  when  the 
discriminant  signals  are  convolved,  using  a  computer,  with  radar  returns  of 
the  target.  When  the  synthesized  discriminant  signals  for  a  preselected 
target  are  convolved  with  the  radar  return  from  a  wrong  target,  the  convolved 
outputs  will  be  significantly  different  from  the  expected  zero  response  or 
single-mode  response,  thus,  the  wrong  target  can  be  discriminated. 

Under  this  research  program,  the  investigation  included  both  theoretical 
analysis  and  experimental  studies  co”,,,''ted  in  a  time-domain  scattering  range. 
Scale  models  of  various  aircrafts  wit;  mplex  geometries  have  been  used  in 
the  experiment.  Over  the  past  few  years,  we  have  largely  accomplished  our 
research  goals.  We  have  proved  the  feasibility  of  our  scheme  and  developed 
it  into  a  workable  and  practical  scheme.  We  feel  that  the  phase  of  basic 
research  has  been  completed  and  it  is  ready  for  us  to  move  into  the  phase  of 
applied  research  and  development. 


1 


2.  RESEARCH  ACHIEVEMENTS 

Since  this  research  program  under  the  sponsorship  of  the  Naval  Air 
Systems  Command  terminates  at  this  point,  it  seems  appropriate  to  review  our 
overall  achievements  of  this  research  program.  They  are  described  in  the 
following  sections: 

2.1  Radar  target  discrimination  by  convolution  of  radar  return  with 
extinction-pulses  and  single-mode  extraction  signals:  This 
section  describes  the  basic  principle  of  our  scheme.  This 
information  has  been  published  in  IEEE  Trans,  on  Antennas  and 
Propagation,  Vol.  AP-34,  No.  7,  pp.  896-904,  July  1986. 

2.2  Radar  target  discrimination  using  the  extinct  ion -pul se  technique: 
This  section  describes  the  extinction-pulse  technique  in  depth.  A 
paper  on  this  subject  has  been  published  in  IEEE  Tran,  on  Antennas 
and  Propagation,  Vol.  AP-33,  No.  9,  pp.  929-  937,  September  1985. 

2.3  Frequency  domain  E-pulse  synthesis  and  target  discrimination:  This 
section  presents  the  frequency  domain  approach  of  our  scheme.  The 
resul  ts  of  this  section  wi  1  1  be  publ  ished  in  the  Apri  1  1987  issue 
of  IEEE  Trans,  on  Antennas  and  Propagation. 

2.4  Extraction  of  the  natural  frequencies  of  a  radar  target  from  a 
measured  response  using  E-pulse  technique:  This  section  describes 
a  new  method  of  natural  frequency  extraction.  This  new  method 
will  be  published  in  a  future  issue  of  IEEE  Trans,  on  Antennas  and 
Propagation. 

2.5  Noise-insensitivity  of  E-pulse  technique  for  target  discrimination: 
This  section  proves  our  scheme  to  be  quite  noise  insensitive.  The 
results  of  this  section  have  been  presented  at  1986  National  Radio 
Science  Meeting,  Philadelphia,  Pennsy  1  van  ia,  June  9-13,  1986. 
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2.6  Determination  of  the  impulse  response  of  a  complex  target 

through  a  deconvolution:  This  section  presents  a  new  method  of 

deconvolution  for  obtaining  the  impulse  response  of  a  target  from 
its  measured  data.  Some  results  of  this  section  have  been 
presented  at  1986  National  Radio  Science  Meeting,  Philadelphia, 
Pennsylvania,  June  9-13,  1986. 
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2.1  RADAR  TARGET  DISCRIMINATION  BY  CONVOLUTION  OF  RADAR  RETURN  WITH 


EXTINCT  ION- PULSES  AND  SINGLE-MODE  EXTRACTION  SIGNALS 


ABSTRACT 


A  new  method  of  radar  target  discrimination  and  identification 
is  presented.  This  new  method  is  based  on  the  natural  frequencies 
of  the  target.  It  consists  of  synthesizing  aspect- independent 
discriminant  signals,  called  Extinction-pulses  (E-pulses)  and 
single-mode  extraction  signals  which,  when  convolved  numerically 
with  the  late-time  transient  response  of  an  expected  target,  lead 
to  zero  or  single-mode  responses.  When  the  synthesized,  discriminant 
signals  for  an  expected  target  are  convolved  with  the  radar  return 
from  a  different  target,  the  resulting  signal  will  be  significantly 
different  from  the  expected  zero  or  single-mode  responses,  thus,  the 
differing  targets  can  be  discriminated.  Theoretical  synthesis  of 
discriminant  signals  from  known  target  natural  frequencies  and 
experimental  synthesis  of  them  for  a  complex  target  from  its 
measured  pulse  response  are  presented.  The  scheme  has  been  tested 
with  measured  responses  of  various  targets  in  the  laboratory. 


U 


INTRODUCTION 


I . 

A  new  radar  discrimination  scheme  has  been  investigated  by  our 
group  over  the  past  few  years  and  this  paper  describes  the  basic 
principle  of  the  scheme  and  reports  some  recent  results.  The  present 
scheme  is  based  on  the  finding  [1-5]  that  for  a  specific  target  there 
exist  aspect- independent  discriminant  signals,  called  the  extinction- 
pulses  (E-pulses*)  and  single-mode  extraction  signals,  which  can  be 
used  to  excite  the  target  to  produce  desirable  late-time  radar  returns 
When  these  aspect-independent  discriminant  signals  with  particular 
waveforms  are  used  to  excite  the  target  they  will  produce'  zero 
response  or  single-mode  responses  in  the  late-time  period.  On  the 
other  hand,  if  the  late-time  response  of  the  target  radar  return, 
which  is  excited  by  a  convenient  radar  pulse  and  is  the  sum  of  target' 
natural  modes,  is  convolved  with  the  discriminant  signals,  the 
convolved  output  will  yield  zero  response  or  single-mode  responses. 

The  former  scheme  is  a  transmitting  scheme  and  the  latter  a  receiving 
scheme.  Since  it  is  difficult  to  physically  synthesize  and  radiate 
the  discriminant  signals  with  particular  waveforms  without  suffering 
distortion,  the  latter  receiving  scheme  is  preferred.  The 
synthesized  signals  with  particular  waveforms  are  now  stored  in  the 
computer;  they  are  convolved  with  the  received  radar  return  inside 
the  computer.  Thus,  we  can  sidestep  the  difficulty  of  synthesizing 
and  radiating  the  discriminalt  signals. 

The  synthesized  discriminant  signals  can  be  used  to  discriminate 
targets  because  when  the  radar  return  of  a  wrong  target  is  convolved 
with  the  synthesized  discriminant  signals  of  the  expected  target, 
the  convolved  outputs  will  be  significantly  different  from  the 


*The  E-pulse  is  similar  to  the  K-pulse  studied  by  other  workers  [5,10] 


expected  zero  response  or  si ng 1 e -mode  responses.  Thus  the  wrong  target  can  be 
discriminated. 

The  complex  natural  resonant  frequencies  of  a  radar  target  are  aspect 
independent  features  of  its  transient  electromagnetic  response.  A  number  of 
researchers  have  recently  attempted  to  discriminate  among  various  targets  by 
extracting  those  natural  frequencies  from  late-time  transient  radar  returns.  Since 
extraction  of  natural  frequencies  from  late-time  target  responses  is  an  inherently 
ill-conditioned  numerical  procedure,  very  large  S-N  ratios  are  required  in  the 
transient  return.  It  has  therefore  been  concluded  that  this  method  for  the  direct 
discrimination  of  differing  targets  is  impractical.  Our  discrimination  scheme 
differs  signi ficantly .  Synthesis  of  the  discriminant  signals  requires  only 
knowledge  of  the  natural  frequencies  of  various  expected  targets.  The  latter 
natural  frequencies  are  measured  in  the  laboratory  where  they  are  extracted  from 
the  late-time  pulse  responses  of  target  scale  models.  The  numerically  ill- 
conditioned  natural  frequency  extraction  procedure  need  therefore  be  applied  only 
to  target  responses  measured  in  a  controlled  S-N  environment.  Synthesized 
discriminant  signals  based  upon  those  laboratory  measurements  are  stored  as  computer 
data  files,  end  subsequently  convolved  numerically  with  actual  transient  target 
radar  returns.  Since  the  latter  convolution  operation  is  numerically  well 
conditioned  (a  smoothing  integral  operator),  the  S-N  requirements  for  the  actual 
radar  return  are  significantly  relaxed. 

Another  observation  made  in  the  course  of  our  study  is  worth  noting.  It  is 
common  thinking  among  many  researchers  that  radar  detection  utilizing  the  late¬ 
time  transient  radar  return  may  not  be  practical  because  it  contains  little  energy; 
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most  energy  is  associated  with  the  early-time  part  of  that  return.  This  thinking 
may  be  true  for  very  low-Q  targets.  Fortunately,  for  most  space  vehicles,  such  as 
rockets  and  aircrafts,  these  targets  are  not  exactly  low-Q  structures.  There  is 
sufficient  energy  contained  in  the  late-time  returns  of  such  targets,  as  can  be 
evidenced  from  our  measured  responses  of  complex  targets  as  discussed  in  this  paper. 


In  Sections  2  and  3,  we  will  demonstrate  how  to  synthesize  theoretically 
and  experimentally  E-0,  E-l,  — /  E-n  pulse  or  the  zero-mode,  the  first-mode, 

- ,  the  nth-mode  excitation  signal.  E-0  pulse  (extinguishing  all  the  redes)  or  the 

zero-mode  extraction  signal  can  be  used  to  produce  a  zero-response  in  the 
convolved  output  while  E-n  pulse  (extinguishing  all  but  the  nth-mode)  or  the  nth- 
mode  extraction  signal  will  produce  a  nth-mode  response  in  the  oonvolved 
output. 

In  Section  4,  examples  are  given  to  show  the  convolution  of  the 

synthesized  discriminant  signals  with  the  measured  radar  returns  of  various 
targets.  The  effectiveness  of  target  discrimination  by  the  present  method  is 
also  discussed. 

Finally,  a  potential  radar  detection  system  based  on  the  present  scheme 
is  suggested  in  Section  5. 

II.  SYNTHESIS  OF  DISCRIMINANT  SIGNALS  3ASED  ON  KNOWN  NATURAL  FREQUENCIES 

The  zero-mode  and  other  single-mode  extraction  signals  or  E-0  pulse  and 
E-n  pulses  of  a  given  target  can  be  synthesized  theoretically  based  on  the 
prior  knowledge  of  natural  frequenices  of  the  target. 

Hie  late-time  radar  return  of  a  target  is  the  sum  of  natural  modes, 
assuming  the  absence  of  noise,  and  it  can  be  expressed  as 
N  at 

h(t)  *  z  a^Q^Je  n  cosfu^t  +  4n(9,$))  (1) 

n  *1 

where  an  is  the  damping  coefficient  and  un  is  the  angular  frequency  of  the  nth 
natural  mode,  and  an(e,<j)  and  $n(9,$)  are  the  amplitude  and  the  phase  angle  of 
the  nth  natural  mode.  It  is  noted  that  a  n  and  u  n  are  independent  of  the 
aspect-angle  (9  and  4),  but  an (9,4 }  and  4(9,4}  are  usually  strong  functions  of 
the  aspect  angle.  N  is  the  total  number  of  natural  modes  which  have  significant 
amplitudes  to  be  considered.  In  practice,  N  can  be  estimated  based  on  the  frequency 
spectrum  of  the  interrogating  radar  pulse. 
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We  aim  to  synthesize  an  extraction  signal  of  duration  Te  which  can  be 
convolved  w:Ja  the  radar  return  h(t)  to  produce  a  single-mode  or  zero-mode 
output  in  the  late-time  period  (t  >  Te) : 


E°(t) 


-/ 

•'A 


E® (t ' )  h(t  -  t')dt',  for  t  >  T* 


where  E°(t)  is  the  convolved  output  signal,  E^t)  is  the  extraction  signal  to 
be  synthesized.  The  substitution  of  (1)  in  (2)  leads  to 
N  V 

E^(t)  *  s  «n(9»*)e  lAn  008  Si*  +  *n(e,*))  +  ^  sin(wnt  +  *n(8,*)}j. 


for  t  >  Tc 


where 


[•  -at* 

6  E^t'J  e  n  cos  unt'  dt' 


re  “V 

Ee(t')  e  sin  wnt'  dt* 


The  coefficients  and  3^  which  determine  the  amplitudes  of  the  natural  modes 

of  the  convolved  output  cure  independent  of  the  aspect  angle  (9  and  i) .  This 

makes  it  possible  to  synthesize  the  aspect- independent  E  (t)  for  producing 

single-mode  or  zero-mode  E°(t) .  It  is  noted  that  An  and  3^  are  numerically 

stable  because  they  are  finite  integrals  over  a  short  period  of  time  T  ,  even 

-at"  e 

though  there  is  a  time  grcwing  factor  of  e  in  them. 

Now  for  example,  if  we  synthesize  Ee(t)  in  such  a  way  that  A^  =  1  and  all 
other  A^  and  go  to  zero,  then  the  output  signal  will  be  a  cosine  first  natural 
mode  as 
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(t)  =  :)  e  cos(^t  +  :>•].  (-,?)),  for  t  >  T0. 

If  we  choose  E^t)  in  such  a  way  that  B3  =  l  and  all  other  ^  and  Bn  go  to 
zero,  then  the  output  signal  will  be  a  sine  third  natural  mode  such  as 

a3t 

E°(t)  =  a30,?)  e  sin(-3t  +  for  t  >  T0. 

If  we  synthesize  E®(t)  in  such  a  way  that  all  ^  and  Bn  vanish,  then  the  zero¬ 
mode  output  is  obtained: 

E°(t)  =*  0,  for  t  >  T0. 

This  Ee'c)  is  the  E-0  pulse  or  the  zero-mode  extraction  signal  because  it 
extinguishes  the  late-time  response  completely. 

The  next  task  is  to  synthesize  E^Ct).  Construct  Ee(t)  with  a  set 
of  basis  function  f^t)  as 
2N 

Ee(t)  *  Z  f^Ct)  (6) 

m=l 

-where  {f  (t) }  can  be  oulse  functions ,  impulse  functions  or  Fourier  cosine 
m 

functions,  etc.  and  {d  ;  are  the  unknown  coefficients  to  be  determined.  We  need  2N  of 

m 

f^t)  because  there  are  2N  of  ^  and  to  be  specified, 

2N  c 

1  Mnm  (7) 

n=l 


Bn 


2N  g 

\  Mnm 
n»l 


where 


_ant’ 

)  e  cos  dt1 


-at* 

)  e  sin  wnt'  dt' 


(8) 


(9) 

(10) 


Equations  (7)  and  (8)  can  be  expressed  in  a  matrix  form  as 
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r  i 

= 

.  s 

®n 

^nm 

n  »  1/2. .. .N 
m  -  1,2. ...2N 

The  coefficient  for  constructing  E0 ( t)  can  then  be  obtained  from 


(11) 


- 

i-  -j 

-1 

«nS 

An 

..  s 

Bn 

(12) 


To  synthesize  £e(t)  for  the  jth-mode  extraction,  we  can  assign  A^  =  1  or  - 

and  let  all  other  A  and  B  go  to  zero,  and  then  calculate  [d_]  from  (12) 

n  n  m 

accordingly.  Using  this  approach,  we  select  an  appropriate  value  for  the 
signal  duration  Tg  which  is  short  and  also  leads  to  a  well  behaved  waveform 
for  Ee(t). 

To  synthesize  the  E-0  pulse  or  the  extraction  signal  for  zero  response, 
all  ^  and  Bn  are  set  to  be  zero.  For  this  case  will  have  non-trivial 


solutions  only  when  the  determinant  of 


det 


C  J 


*4 


vanishes.  That  is 


(13) 


This  condition  can  be  met  because  all  the  elements  of  Mc  and  Ms  are  functions 

nm  nm 

of  the  signal  duration  T  ,  and  it  is  possible  to  numerically  search  for  the 

optimum  value  of  Tg  which  makes  (13)  valid.  Usually  we  use  the  Newton  method 

for  this  purpose.  Once  the  optimum  Tg  is  determined,  [dm]  can  be  easily 

determined  from  a  set  of  homogeneous  equations  generated  from  (11)  by  setting 

A  and  B  to  be  zero, 
n  n 

Single-mode  extraction  signals  or  E-n  pulses  can  also  be  synthesized  in  a 
S':"ilar  way  as  that  for  the  E-0  pulse.  For  example,  if  we  aim  to  synthesize  the 
cosine  first  mode  extraction  signal,  we  drop  the  first  equation  involving 
A^  in  (11).  We  then  have  a  new  set  of  (2N-1)  simultaenous  equations  as 
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(14) 


~A2' 

..  c 

/ 

“  — 

Mnm 

= 

B1 

„  s 

Mnm 

_®N_ 

_ 

n  =  2,3 - N  for 

n  =*  1,2. . .  .N  for  M_f 
m  -  1,2 _ 2N-1 


where 

7  is  a  (2N-1)X(2N-1)  matrix  and  it  is  the  original 

,Mnm  . 

• 

.«ri. 

matrix  with  its  first  row  and  its  last  column  removed.  We  now  force  A2 , 
Ajj,  B^,  ....  Bjj  to  be  zero.  This  requires  that 


det 


(15) 


Equation  (15)  can  be  solved  numerically  to  determine  the  optimum  Tg 
for  the  cosine  first  mode  extraction  signal.  After  that  [d  ]  can  be 
determined  from  a  corresponding  set  of  homogeneous  equations  from  (14) . 

It  has  been  found  that  the  signal  duration  Tg  for  various 
discriminant  signals  is  not  unique  but  there  exists  an  optimum  value  of 
Te  for  which  the  waveform  of  the  discriminant  signal  is  best  behaved 
and  it  possesses  a  maximal  sensitivity  in  discriminating  the  target  ’6] . 
This  optimum  value  of  Tq  is  numerically  found  to  be  slightly  longer 
than  twice  the  transit  time  of  the  target  or  2  L/C  where  L  is  the 
longest  dimension  of  the  target  and  C  is  the  speed  of  light. 


III.  SYNTHESIS  OF  DISCRIMINANT  SIGNALS  BASED  ON  EXPERIMENTAL  PULSE 

RESPONSE 

In  the  preceding  section,  discriminant  signals  for  zero  response 
or  single-mode  responses  were  synthesized  based  on  the  prior  knowledge 
of  target's  natural  frequencies.  In  the  case  of  complex  targets,  this 
information  is  difficult  to  obtain,  and  the  synthesis  of  the  discriminant 
signals  will  be  based  on  an  experimental  measurement  of  the  pulse 
response  or  the  late-time  response  of  the  scale  model  of  the  target 
combined  with  some  theoretical  techniques  described  below. 


A.  Application  of  Continuation  Method  to  a  Measured  Pulse  Response 

The  first  method  we  have  employed  to  synthesize  the  discriminant  signals  for 
a  target  is  to  apply  a  recently  developed,  continuation  method  [7]  to  extract 
major  natural  frequencies  of  the  target  frcm  its  measured  pulse  response. 

This  method  is  to  regularize  the  ill-conditioned  problem.  After  the  natural 
frequencies  are  determined,  the  theoretical  technique  of  the  preceding  section 
is  used  to  synthesize  the  discriminant  signals.  The  following  steps  are  used  in 
the  process: 

(1)  Measure  the  pulse  response  of  the  scale  model  of  the  target  at  various 
aspect  angles. 

(2)  Use  the  Feist  Fourier  Transform  to  obtain  approximate  values  of  natural 
frequencies  frcm  the  measured  pulse  response. 

(3)  Employ  the  continuation  method  and  the  natural  frequencies  obtained  frcm  FFT 
as  the  initial  guesses  to  calculate  accurate  values  of  natural  frequencies  of 
the  target.  At  each  step  of  the  algorithm ,  the  condition  number  of  the 
regularized  problem  is  checked. 

(4)  The  theoretical  technique  of  the  preceding  section  is  then  applied  to  synthesize 
the  discriminant  signals  for  the  target. 

B.  Application  of  Fast  Prony's  Method  to  a  Measured  Pulse  Response 

The  second  method  we  have  used  to  synthesize  the  discriminant  signals  for  a  target 
is  mainly  based  on  the  fast  Prony's  method  [3,8]  and  a  measured  pulse  response.  This 
method  in  discrete-time  basis  consists  of  the  following  steps: 

(1)  The  fast  Prony's  method  is  applied  to  the  sampled  data  of  the  measured  pulse 
response  of  the  target.  The  coefficients  of  Prony  polynomial  are  obtained. 

The  time-ordered  sequence  of  these  coefficients  can  be  shown  to  became 
the  E-0  pulse  [E-0]  of  the  target  [8] ,  which  eliminates  the  natural  modes 
of  the  target  and  other  noise. 
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(2)  Find  coots  of  the  Z-transform  of  the  sequence  which  respresents  [2-0]  . 

From  those  roots  in  the  proper  locations  of  the  complex  plane,  the 
complex  natural  frequencies  [Sm]  of  the  target  can  be  calculated. 

(3)  At  this  step,  E-0  pulse  [E-0]  which  extinguishes  N  natural  modes  of  the  target 
can  be  constructed  by  convolving  N  couplets  (1,  -z  )  for  m  =  1,2,...N  wh.^ore 

=  exp(SmlT)  and  IT  is  the  sampling  interval: 
[E-0]«(1,-z1)*(1,-z2)*....*(1,-zn) 

(4)  The  j th-mode  extraction  signed  CEje]  can  be  obtained  by  removing  the 
couplet  (l,-Zj)  from  [2-0]  as 

[Ejl  ■  (l,-z1)*....*(l,-zi)*(l,-Z|c)*....*(l,-zN) 

4.  CONVCLUTICN  OF  SYNTHESIZED  DISCRIMINANT  SIGNALS  WITH  EXPERIMENTAL  RADAR 
RETURNS 

To  demonstrate  the  applicability  of  the  present  target  discrimination 
method  in  practice,  we  have  convolved  the  measured  radar  returns  of  various 
targets  with  the  corresponding  synthesized  signals  for  extracting  zero-response 
and  single-mode  responses.  We  have  tested  many  targets  and  produced  many 
interesting  results  but  only  a  few  examples  will  be  given  here  for  brevity. 

Experimental  pulse  responses  of  various  radar  targets  used  in  the 
present  study  were  either  measured  recently  in  the  time-domain  scattering  ran^e 
at  Michigan  State  University  or  have  been  published  previously  [3].  The  MSU 
scattering  range  uses  a  biconical  transmitting  antenna  (2.5  m  long,  16°  cone 
angle)  and  as  a  receiving  probe  it  uses  a  short  monopole  field  probe  on  the 

ground  plane  (6  m  x  5  m),  a  current  probe  or  a  charge  probe  on  the  target  surface. 
The  transmitting  antenna  is  excited  either  by  a  mercury-switched  nanosecond 

pulser  (Tektronix-109)  ,  which  produces  pulses  of  about  lOOpS  risetime  and 
variable  duration  1  ns  <  tp  <  1  us  with  amplitude  as  great  as  500  volts  at  a 
1  kHz  rate,  or  by  another  picosecond  pulser  (Picosecond  Lab-IOOOB)  which  produces 
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a  narrower  but  a  lower  amolitude  oulse.  T'ne  signal  from  the  receiving  probe  is 
measured  by  a  sampling  scope  and  the  signal  processing,  including  averaging  and 
clutter  subtraction,  is  performed  by  a  computer  controlled  system. 

Figure  1  shows  the  measured  impulse  response  of  a  thin  cylinder  target, 
having  length  t  *  12.5  in  »  31.75  cm  and  length-to-radius  ratio  1/ a  *  400, 
illuminated  obliquely  by  a  nanosecond  plane-wave  pulse  at  8  3  60°  aspect  and 
the  synthesized  E-0  pulse.  This  E-0  pulse  was  synthesized  in  terms  of  pulse 
functions  with  the  first  three  natural  modes ,  indicated  in  the  figuro,  extracted 
by  the  continuation  method  as  described  in  Section  3. A.  The  band- limited  frequency 
spectrum  of  tha  1-nS  incident  pulse  precluded  the  indentification  of  higher  natural 
frequencies.  Figure  2  shows  the  convolved  result  of  the  E-0  pulse  of  Fig.  1 
with  the  pulse  response  of  the  correct  cylinder  along  with  the  similar 
convolved  results  with  the  pulse  responses  of  cylinders  5%  and  25%  longer. 

The  correct  cylinder  yields  a  low-amplitude  late-time  response  which 
approximates  the  expected  zero  response?  failure  to  achieve  a  perfect  zero 
response  is  the  result  of  imperfect  experimental  data  having  content  other 
than  the  assumed  sum  of  pure  natural  modes  (the  latter  may  be  artifacts 
introduced  by  the  measurement  system).  Longer  targets  lead  to  responses  which 
can  be  discriminated  from  an  expected  zero  response.  Discrimination  of  the  5% 
longer  cylinder  is  marginal  while  that  of  the  25%  longer  target  is  clear. 

Next  example  was  given  to  show  the  convolution  of  single-mode  extraction 
signals  with  the  measured  radar  return  of  a  target.  Figure  3  shows  the 
results  of  the  convolution  of  the  measured  pulse  response  of  a  wire  target 
with  18.6"  length  and  .0465"  radius  with  the  synthesized  signal  for  the  third 
mode  extraction.  Figure  3(a)  shows  the  synthesized  signal  for  extracting  the 
cosine  third  mode,  constructed  in  terms  of  impulse  basis  functions  using  the 
fast  Prony's  method  described  in  Section  33,  and  the  convolved  output.  It  is 
observed  that  a  cosine  third  mode  is  produced  in  the  late-time  period  of  the 


convolved  output.  Figure  3(b)  shows  the  synthesized  signal  for  extracting  the 
sine  third  mode  (in  terms  of  impulse  basis  functions)  and  convolved  output.  A 
sine  third  mode  is  produced  in  the  late-time  period.  To  demonstrate  that 
these  late-time  responses  are  indeed  that  of  the  third  mode,  a  complex 
convolved  output  is  constructed  as 
C(t)  *  A(t)  -  j  B (t) 

where  A(t)  is  the  convolved  output  of  Fig.  3(a)  and  B(t)  is  the  convolved 
output  of  Fig.  3  (b).  If  the  late-time  response  of  A(t)  is  a  pure  cosine 
third  mode,  it  can  be  expressed  as 

a3t 

[AttJl^g  3  a3  e  costa^t  +  $3). 

Similarly,  if  the  late-time  response  of  B(t)  is  a  pure  sine  third  mode,  it  can 
be  expressed  as 

°3t 

*  a3  e  sin<w3t+$3). 

Thus,  the  late-time  response  of  the  complex  convolved  output  should  be 

a  3 1  -jto^t  +  $  3 ) 

[c(t)1late  =  a3  e  e 

The  logarithm  of  the  late-time  response  of  C(t)  yields 
ln[C(t)llate  »  In  a3  +  <J3t  -  jfa3t  +  $3) 
and 

real  part  of  ln[C(t)Jlate  =  <73t  +  In  a3 
imaginary  part  of  lnlCtt)]^^  ■  -u3t  -  $3 
Figure  3(c)  shows  the  real  part  of  ln[C(t)J  in  comparison  with  the  line  of 
(a3t).  It  is  observed  that  the  late-time  response  of  ReQnICCt)])  is  nearly  in 
parallel  with  the  line  of  (a^t)  for  the  late-time  period  of  t/(L/C)  >  3  where 
t/(I JO  is  the  normalized  time  with  L  as  the  wire  length  and  C  as  the  speed  of 
light.  Figure  3(d)  shows  the  imaginary  part  of  ln(C(t)J  in  comparison  with 
the  line  of  (-«•>  t).  It  is  again  observed  that  the  late-time  response  of 
ImUnlC(t)]}  is  in  parallel  with  the  line  of  (-u>3t)  for  the  late-time  period 
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of  t/(L/C)  >  3.  Figures  3(c)  and  3(d)  give  a  positive  indication  that  the 
late-time  response  of  the  convolved  outputs  of  Figs.  3(a)  and  3(b)  contain 
only  (or  predominantly)  the  third  mode  of  the  target. 

An  example  is  also  aiven  to  show  the  capability  of  tarqet  discrimination 
provided  by  the  present  method.  Figure  4  shows  the  results  of  the  convolution 
of  the  radar  return  of  a  wrong  target,  a  16.74"  wire,  with  the  synthesized  siqnal 
for  extracting  the  third  mode  of  the  right  target,  a  18.6"  wire.  In  Figs.  4(a) 
and  4(b),  the  convolved  outputs  are  significantly  different  from  the  expected 
third  mode  of  the  right  tarqet.  Furthermore,  the  late  time  response  of  Re  {ln[C(t)]> 
deviates  greatly  from  the  line  of  (c^t)  and  the  late-time  response  of  Im  (ln[C(t)]} 
also  differs  significantly  from  the  line  of  (-u^t)  as  shown  in  Figs.  4(c)  and 
4(d).  The  results  given  in  Fig.  4  give  a  positive  discrimination  of  the  wrong 
target,  implying  that  the  radar  return  does  not  belong  to  the  right  target  of 
18.6"  wire. 

We  have  also  applied  the  present  method  to  discriminate  between  complex 
radar  targets.  Two  tarqets,  experimental  models  of  the  F-18  plane  and  the  7 07 
plane  which  have  similar  sizes  but  different  geometries,  were  used  for  this 
purpose.  In  the  experiment,  the  half  of  each  model  was  placed  on  the  ground  plane 
(image  effect),  and  was  illuminated  by  a  pulsed  electric  field  polarized 
perpendicularly  to  the  ground  plane.  Results  on  these  two  targets  are  depicted 
in  the  following  figures. 

Figure  5(a)  shows  the  convolution  of  the  E-0  pulse  for  the  F-18  plane  with 
the  measured  radar  response  of  the  F-18  Diane.  The  convolved  outputs  shows  a 
strong  early-time  response  followed  by  a  "extinguished"  (almost  zero)  late-time 
response  as  expected.  Figure  5(b)  shows  the  convolved  output  of  the  E-0  pulse 
for  the  707  plane  with  the  measured  radar  response  of  the  F-18  plane.  This 
convolved  output  shows  a  significant,  "unextinguished"  late-time  response  imnlvin^ 
that  the  707  E-0  pulse  was  convolved  with  the  radar  response  of  a  wrong  target 
other  than  the  707  plane. 
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Figure  6  shows  the  convolved  results  of  the  measured  radar  resoonse  of  the 
707  plane  with  the  fourth-mode  extraction  signal  for  the  707  plane.  In  this  figure, 
the  extracted  angular  frequency  line  is  almost  parallel  to  the  ~4t  line  (of  the 
707)  and  the  extracted  damping  coefficient  line  closely  parallel  to  the  r4t  line 
(of  the  707)  in  the  late-time  period.  This  implies  that  the  radar  response 
belongs  to  the  right  target  of  the  707  plane. 

Figure  7  shows  the  convolved  results  of  the  measured  radar  response  of 
the  F-18  plane  with  the  fourth-mode  extraction  signal  for  the  707  plane.  The 
extracted  angular  frequency  line  deviates  from  the  (co^t)  line  of  the  707  plane 
and  the  extracted  damping  coefficient  line  also  differs  from  the  (c.t)  line  of 
the  707  plane.  These  results  imply  that  the  radar  response  belongs  to  a  wronq 
target  other  than  the  707  plane. 

It  is  noted  that  the  synthesized  E-0  pulses  and  single-^ode  extraction 
signals  for  these  targets  are  available  elsewhere  [9], 

5.  A  POTENTIAL  RADAR  DETECTION  SYSTEM  BASED  CN  THE  E8ESENT  OONVCLOTION 

METHOD 

Figure  8  depicts  a  potential  radar  detection  system  based  on  the  present 
scheme  of  convolving  the  target  radar  return  with  the  synthesized  discriminant 
signals.  The  system  consists  of  a  network  of  computers  and  each  of  them  is 
assigned  to  store  the  synthesized  signals  for  extracting  various  single-mode  or 
zero-mode  response  for  a  particular  friendly  target.  All  the  relevant 
friendly  targets  are  assumed  to  be  covered  in  the  network  of  computers.  When 
an  approaching  target  is  illuminated  by  an  interrogating  radar  signal,  the 
radar  return  is  divided  and  fed  to  each  computer  after  amplification  and 
signal  processing.  Inside  each  computer  the  stored  discriminant  signals  are 
convolved  with  the  radar  return.  In  principle,  only  one  of  the  computers  will 
produce  various  single-mode  and  zero-mode  outputs  in  the  late-time  period;  the 
rest  of  the  computers  should  produce  irregular  outputs.  The  computer 
producing  the  single-mode  and  zero-mode  outputs  will  then  be  identified  with 
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the  target.  If  none  of  the  computers  produces  single-mode  and  zero  outputs 
in  the  late-time  period,  the  approachinq  tarqets  will  not  be  a  friendly  tarqet. 

6.  Discussion 

The  present  radar  discrimination  scheme  is  not  a  conventional  correlation 
method  because  the  measured  radar  return  is  convolved  with  synthesized  discriminant 
signals  which  possess  particular  waveforms.  The  present  scheme  is  basically 
different  from  the  matched  filterinq  scheme.  A  fundamental  difference  between 
the  E-pulse  scheme  and  the  matched  filtering  scheme  is  as  follows:  The  E-pulse 
scheme  is  to  nullify  a  selected  natural  mode  set,  and  does  not  deal  with 
maximizing  a  specific  modal  content.  For  example,  the  E-0  pulse  eliminates  all 
the  modes  in  the  late  time.  The  single-mode  extraction  signal  eliminates  all 
but  one  mode.  We  do  not  consider  maximization  of  responses  in  the  presence  of 
noise  at  specific  times  as  in  the  matched  filtering  scheme.  A  further  difference 
between  these  two  schemes  is  the  following:  The  natural  modes  are  not,  in  general, 
orthogonal.  Maximizing  the  response  of  a  particular  mode  will  not  nullify  the 
other  modes.  Also  note  that  the  E-pulse  scheme  is  aspect-independent,  whereas 
the  matched  filtering  is  not. 

We  have  not  attempted  to  compare  the  performance  of  the  E-pulse  scheme  with 
that  of  a  matched  filtering  scheme. 
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Synthesis  of  Thin-Cyl inder  E-0  Pulse  from  Experimentally 
Measured  Pulse  Response;  Dominant  Natural  Frequencies 
Extracted  by  the  Continuation  Method. 
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Fig.  3.  Convolution  of  the  impulse  response  of  a  18.6"  (length)  wire 
with  the  synthesized  signal  for  the  third  mode  excitation 
of  the  18.6"  wi re. 


23 


4.  Convolution  of  the  impulse  response  of  a  wrong  target,  a  16.74" 
wire,  with  the  synthesized  signal  for  exciting  the  third  mode 
of  the  right  target,  a  18.6"  wire. 
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gular  frequency  and  damping  coefficient  extracted  from  the  complex  con 
de  excitation  signal  for  the  707  target  and  the  F-18  measured  response 
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Fig.  8.  A  proposed  target  identification  and  discrimination  system. 


2.2  RADAR  TARGET  DISCRIMINATION  USING  THE  EXTINCTION- PULSE  TECHNIQUE 


ABSTRACT 

An  aspect  independent  radar  target  discrimination  scheme  based  on  the 
natural  frequencies  of  the  target  is  considered.  An  extinction-pul se  waveform 
upon  excitation  of  a  particular  conducting  target  results  in  the  elimination 
of  specified  natural  modal  content  of  the  scattered  field.  Excitation  of  a 
dissimilar  target  produces  a  noticably  different  late-time  response. 
Construction  of  appropriate  extinction-pulse  waveforms  is  discussed,  as  well 
as  the  effects  of  random  noise  on  their  application  to  thin  cylinder  targets. 
Also  presented  is  experimental  verification  of  this  discrimination  concept 
using  simplified  aircraft  models. 
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Introduction 


Radar  target  identification  methods  using  the  time-domain  response  of  a 
target  to  a  transient  incident  waveform  have  gained  considerable  interest 
recentl y  [ 1-4] .  One  of  the  most  intriguing  schemes  involves  the  so-called 
"Kill-pulse"  technique  as  first  described  by  Kennaugh  [5].  A  Kill-pulse  (K- 
pulse)  is  an  excitation  waveform  synthesized  in  such  a  way  as  to  minimize  a 
transient  scattered  field  response.  Target  discrimination  results  from  the 
unique  correspondence  of  a  K-pulse  to  a  particular  target;  excitation  of  a 
dissimilar  target  yields  a  "larger"  response. 

This  paper  describes  a  related  "extinct ion- pul se"  (E-pulse)  concept, 
based  on  the  natural  resonance  of  a  conducting  radar  target  via  the 
singularity  expansion  method  (SEM)  [6].  The  time  domain  electric  field 
scattered  by  the  target  is  divided  into  an  early-time,  forced  response  period 
when  the  excitation  waveform  is  traversing  the  target,  and  a  late-time,  free 
oscillation  period  that  exists  after  the  excitation  waveform  has  passed  [7,8]. 
The  early-time  response  is  not  utilized  due  to  its  complicated  nature.  The 
late-time  response  can  be  decomposed  into  a  sum  of  damped  sinusoids 
oscillating  at  frequencies  determined  entirely  by  the  geometry  of  the  target. 
An  E-pulse  is  then  viewed  as  a  transient,  finite  duration  waveform  which 
anihilates  the  contribution  of  a  select  number  of  these  natural  resonances  to 
the  late-time  response.  A  related  target  identification  scheme  based  on 
natural  target  resonances  has  been  examined  by  Chen  [9]. 

Target  discrimination  using  this  SEM  viewpoint  is  easily  visualized. 

Each  target  can  be  described  by  a  set  of  natural  frequencies.  An  E-pulse 
designed  to  annul  certain  natural  resonances  of  one  target  will  excite  those 
of  a  target  with  different  natural  frequencies,  resulting  in  a  different 
scattered  field.  Also  made  apparent  is  the  aspect-independent  nature  of  the 
E-pulse.  Since  the  values  of  the  target  resonance  frequencies  are  independent 
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of  the  excitation  waveform,  the  E-pulse  will  eliminate  the  desired  natural 
modal  content  of  the  late-time  scattered  field  regardless  of  the  orientation 
of  the  target  with  respect  to'  the  transmitting  and  receiving  antennas. 

It  is  important  to  note  that  the  E-pulse  waveform  need  not  be  transmitted 
to  employ  this  concept.  It  is  assumed  that  an  excitation  waveform  with  finite 
usable  bandwidth  will  be  used  to  excite  the  target,  resulting  in  a  measured 
scattered  field  with  the  desired  (finite)  modal  content.  The  E-pulse  can  then 
be  convol  ved  numerical  ly  with  the  measured  target  response,  yielding  results 
analogous  to  E-pulse  transmi ssion.  If  the  maximum  modal  content  of  the  target 
scattered  field  car,  be  estimated  from  the  frequency  content  of  the  excitation 
pulse,  then  the  E-pulse  waveform  can  be  constructed  to  yield  a  null  late-time 
convolved  response. 

After  an  initial  presentation  of  the  SEM  representation  of  the 
backscattered  field  excited  by  a  transient  incident  wave  and  calculation  of 
the  corresponding  impulse  response,  two  types  of  E-pulses,  forced  and  natural, 
will  be  discussed.  The  results  are  then  specialized  to  a  thin  cylinder 
target,  which  has  an  impulse  response  amenable  to  analytic  calculation. 

Target  discrimination  using  the  natural  thin  cylinder  E-pulse,  as  well  as  the 
effects  of  random  noise  are  also  investigated.  Lastly,  experimental 
verification  of  the  E-pulse  concept  is  presented. 

II.  Backscattered  Field  Excited  by  Transient  Incident  Wave 

A  perfectly  conducting  radar  target  is  illuminated  by  a  plane 
electromagnetic  wave  as  shown  in  Figure  1.  The  electric  field  associated  with 
this  transient  wave  can  be  written  in  the  Laplace  transform  domain  as 
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(1) 


—  k  *F 

V{r,s)  =  ;e(s)e  c 

where  r  is  a  position  vector  in  a  coordinate  system  local  to  the  target,  ;  is  a 
constant  unit  vector  specifying  the  polarization  of  the  wave,  k  is  a  unit 
vector  in  the  direction  of  propagation,  and  e(t)  represents  the  time  dependence 
of  the  incident  field.  The  current  £  induced  on  the  surface  of  the  target  is 
given  by  the  solution  to  the  transform  domain  E-field  integral  equation 


jC 


2  „ 

?*iT(r', s)(t*7)  -  t*k(r',s) 

c^ 


-1* 
,  c 
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ds  *  =  -Se0t ‘E1 (r,s) 


(2) 


where  t  is  a  unit  vector  tangent  to  the  surface  of  the  target  at  the  point  1* 
on  S,  and  R  s|  r-r’l.  The  surface  current  can  also  be  expanded  using  the  SEM 


representation 


it(f,s)  =  I  a2Ka(r)(s-s3)“1  +  W(r,$) 


a-1 


(3) 


where  W(r,s)  represents  any  entire  function  contribution  to  the  current,  and 
it  is  assumed  that  only  a  finite  number  of  natural  modes  are  substantially 
excited  by  the  incident  field.  It  is  further  assumed  that  the  only 
singularities  of  K(r,s)  in  the  finite  complex  s-plane  are  simple  poles  at 
natural  frequencies  s^  =  Then,  K  (r)  represents  the  current 

distribution  of  the  tfh  natural  mode,  and  a^  is  the  "class-1"  coupling 
coefficient  given  by  Baum  [10 ]. 

The  far-zone  backscattered  electric  field  can  be  computed  by  integrating 

the  current  distribution  over  the  suface  of  the  target.  The  ;  component  of 

backscattered  field  can  then  be  written  as 

s 

— r 
„  c 
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e(s)H(s,k,; ) 
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where  r  =  ,  r  ,  and  H(s,k,;)  is  the  aspect  dependent  transfer  function  of  the 
target.  Using  (3)  to  represent  the  surface  current,  the  transfer  function  can 
be  inverse  transformed  to  determine  the  backscatter  impulse  response  of  the 
target . 

Because  of  the  entire  function  contribution  to  the  current,  the  impulse 
response  will  exhibit  two  distinct  regions.  The  early-time,  forced  component 
represents  the  backscattered  field  excited  by  currents  during  the  time  when 
the  impulse  is  traversing  the  target;  it  has  a  duration  equal  to  twice  the 
one-way  maximal  transit  time  of  the  target,  T.  The  late-time  free  oscillation 
component  is  composed  purely  of  a  sum  of  constant  amplitude  natural  modes  and 
exists  for  all  time  t  >  2T  as 

N  ,  -  o  t 

h(t,k,0  =  z  an(k,;)e  cos(^nt  +  <f>(k,;))  t  >  2T  (5) 

n=l 

where  it  has  been  assumed  that  the  entire  function  makes  no  contribution  to 
the  late-time  component  [6].  Thus,  the  ;  component  of  the  far-zone 
backscattered  field  is  given  simply  by  the  convolution  of  the  time  domain 
incident  field  and  the  impulse  response,  and  is  also  composed  of  forced  and 
freely  oscillating  portions. 


III.  Incident  E-pulse  Waveform  Synthesis 

To  synthesize  an  E-pulse  for  a  particular  target,  the  convolutional 
representation  of  the  backscattered  field  is  written  in  integral  form  using 
the  impulse  response  of  (5)  as 


E?(t 


e (t ' )h (t-t ' ,k,;)dt' 


e  ..  ^  ^  an( t-t ' ) 


[  Te  N 
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This  response  is  valid  for  the  late-time  portion  of  the  scattered  field,  t  > 
TL  =  Tg  +  2T,  where  Te  is  the  duration  of  e(t).  The  excitation  waveform 
becomes  an  E-pulse  when  the  scattered  field  is  forced  to  vanish  identically 
the  late-time.  Rewriting  equation  (6)  and  employing  this  condition  yields  a 
defining  equation  for  the  E-pulse 


5n^ 


s  an0qq)e  [An(Te)cos(^nt  +  ;n(k,;))  +  Bn (Te)s1  n (--int  +  >n(k,;))]  »  0 
n  =1 

t  >  \  -  Te  +  2T  (7) 

where  the  coefficients  An(Te)  and  Sn(Te)  are  given  by 
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(8) 


The  linear  independence  of  the  damped  sinusoids  in  equation  (7)  requires 


An(Te)  =  Bn(Te)  =  0  for  all  1  <  n  <  H. 


It  is  important  to  note  that  An(Te)  and  Bn(Tg)  are  independent  of  the 
aspect  parameters  k  and  verifying  the  aspect  independence  of  the  E-pulse. 
This  is  a  direct  consequence  of  the  separability  of  the  terms  of  the  impulse 
response . 

A  physical  interpretation  of  the  E-pulse  can  be  facilitated  by 
decomposing  the  excitation  waveform  as  shown  in  Figure  2  as 

e(t)  =  ef(t)  +  ee(t)  (9) 

where  e^(t)  is  an  excitatory  component,  nonvanishing  during  0  <t  <Tf,  the 
response  to  which  is  subsequently  extinguished  by  ee(t)  which  follows  during 
Tf  <  t  <  Tg.  Substituting  equation  (9)  into  equation  (8)  and 
using  An(Te)  =  8n(TQ)  *  0  yields 


f^'-v  Icos^f; 

er(t')e  j  , 

(sin^t1. 


dt' 


(10) 
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The  extinction  component  of  the  E-pulse  necessary  to  eradicate  the  response 
due  to  a  preselected  excitatory  component  can  be  constructed  as  an  expansion 
over  an  appropriately  chosen  set  of  linearly  independent  basis  functions  as 

ee^)  =  z\  dD 

m=1 


Equation  (10)  then  becomes 


2N 
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1  <  i  <  N 


(12) 


where 


This  can  be  written  using  matrix  notation  as 


(13a) 

(13b) 


Solving  this  equation  for  Cj_,  C2N  determines  the  extinction  component  via 
equation  (11)  and  thus  the  E-pulse. 

It  is  convenient  at  this  point  to  identify  two  fundamental  types  of  E- 
pulses.  When  T^  >  0  the  forcing  vector  on  the  right  hand  side  of  equation 
(14)  is  non-zero  and  a  solution  for  ee(t)  exists  for  almost  all  choices  of  Tg. 
This  type  of  E-pulse  has  a  non-zero  excitatory  component  and  is  termed  a 
"forced"  E-pulse.  In  contrast,  when  Tf  *  0  the  forcing  vector  vanishes  and 
solutions  for  ee(t)  exist  only  when  the  determinant  of  the  coefficient  matrix 
vanishes,  i.e.,  when  detCMfrg)]*  0.  These  solutions  correspond  to  discrete 
eigenvalues  for  the  E-pulse  duration  T  ,  which  are  determined  by  rooting  the 
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determi nantal  characteristic  equation.  Since  there  is  no  excitatory 
component,  this  type  of  E-pulse  is  viewed  as  extinguishing  its  own  excited 
field  and  is  called  a  "natural"  E-pulse. 

IV.  Thin  Cylinder  E-pulse  Analysis 

A  theoretical  analysis  of  a  thin  wire  target  has  been  undertaken  by 
various  authors  [11,12].  Target  natural  frequencies  are  determined  from  the 
homogeneous  solutions  to  the  integral  equation  (2).  The  geometry  of  the 
target  and  its  orientation  with  respect  to  the  excitation  field  are  given  in 
Figure  3  along  with  the  first  ten  natural  frequencies.  The  frequencies  are 
normalized  by  *c/L  where  L  is  the  length  of  the  wire  and  c  is  the  speed  of 
light,  and  correspond  to  a  wire  of  radius  given  by  L/a  =  200. 

The  thin  cylinder  impulse  response  can  be  calculated  by  inverting 
equation  (4)  [12]  and  becomes  a  pure  sum  of  natural  modes  in  the  late'time. 
Figure  4  shows  the  impulse  responses  of  thin  cylinders  oriented  at  3  =  30°  and 
■3  =  60°,  generated  by  using  the  first  ten  modes  of  the  target.  Note  the 
distinct  early  and  late-time  regions.  An  E-pulse  to  extinguish  the  first  ten 
modes  of  the  target  is  created  by  solving  equation  (14)  using  the 
corresponding  frequencies. 

Figure  5  shows  a  natural  E-pulse  synthesized  to  kill  the  first  ten 
natural  modes  of  the  thin  wire  target,  using  a  pulse  function  basis  set.  The 
duration  of  the  waveform,  Te  =  2.0408  l/c,  corresponds  to  the  first  root  of 
the  resulting  determinantal  equation.  Superimposed  with  this  is  Kennaugh’s 
original  K-pulse.  The  similarlity  is  striking,  with  the  major  difference 
being  the  finite  duration  of  the  E-pulse.  Also  shown  in  Figure  5  is  a  forced 
pulse  function  E-pulse  constructed  to  extinguish  the  first  ten  modes  of  the 
target.  The  duration  has  been  chosen  as  Tfi  =  2.3  and  the  excitation  component 
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has  been  chosen  as  a  pulse  function  of  width  equal  to  that  of  the  basis 
functions . 

Numerical  verification  of  the  thin  wire  E-pulse  is  given  in  Figure  6. 

The  natural  E-pulse  waveform  of  Figure  5  is  convolved  with  the  30°  and  60° 
impulse  responses  of  Figure  4  and  the  resulting  backscattered  field 
representations  are  observed  to  be  zero  in  the  late-time.  Note  also  the 
expected  non-zero  early-zime  response.  This  portion  is  useful  since  it 
provides  a  comparative  benchmark  assessing  the  quality  of  the  annulled 
component  of  the  response;  this  is  important  when  imperfect  "excti notion"  is 
evident  (due  to  noise,  errors  in  the  natural  frequencies,  etc.). 

The  question  of  E-pulse  waveform  uniqueness  as  the  number  of  natural 
frequencies  extinguished,  N,  is  taken  to  infinity  has  not  been  resolved.  For 
the  thin  cylinder  target  there  appears  to  be  such  a  unique  waveform.  Evidence 
is  provided  by  Figure  7  which  shows  natural  E-pulses  of  minimum  duration 
designed  to  extinguish  various  numbers  of  natural  modes,  using  both  Fourier 
cosine  and  pulse  function  basis  sets.  It  is  apparent  that  these  waveforms  are 
nearly  identical  and  they  appear  to  be  converging  to  a  particular  shape. 

The  most  critical  parameter  necessary  for  E-pulse  uniqueness  is  the 
finite  E-pulse  duration,  Tg.  This  must  converge  to  a  distinct  value  as  N—»  . 
For  the  pulse  function  basis  set  it  can  be  shown  that  the  natural  E-pulse 
durations  are  given  "by 

T  *  2*N  P  =  1,2, ...  1  <  z  <  N  (15) 

i 

where is  the  imaginary  part  of  the  nth  natural  frequency.  If  this  is 
written  for  the  thin  cylinder  as 

-i  »  -  5(0]  l  (16) 
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where  5(0/ i  decreases  ass.ymptotical  ly  to  zero  with  increasing  then  the 
minimum  E-pulse  duration  converges  to 


1  i  m  2-N  L 
( 'e'min  '  -[ N-5 ( n) ]  c 

-2L- 

c 


(17) 


V.  Target  Discrimination  with  E-pulse  Waveforms 

Discrimination  between  different  thin  cylinder  targets  is  demonstrated 
by  convolving  the  natural  E-pulse  of  Figure  5,  which  has  been  constructed  to 
extinguish  the  first  ten  modes  of  a  target  of  length  L,  with  the  impulse 
response  of  the  expected  target  and  a  target  5%  longer.  The  result  is  shown 
in  Figure  8.  The  late-time  response  of  the  expected  target  has  been 
successfully  annulled,  while  the  response  of  the  differing  target  is  non-zero 
over  the  same  period.  The  difference  in  target  natural  frequencies  provides 
the  basis  for  discrimination  based  on  the  comparison  of  adequatel y  di ssimi 1  ar 
late-time  responses  of  differing  targets. 

Sensitivity  of  E-pulse  performance  to  the  presence  of  uncorrelated 
random  noise  is  investigated  by  perturbing  each  point  of  the  thin  cylinder 
impulse  response  of  Figure  4  by  a  random  amount  not  exceeding  10%  of  the 
maximum  waveform  amplitude.  The  result  is  shown  in  Figure  9.  An  attempt  is 
then  made  to  extinguish  this  noisy  response  by  convolving  it  with  the  natural 
E-pulse  of  Figure  5.  As  expected,  the  convolution,  shown  in  Figure  10,  does 
not  exhibit  a  null  late-time  response,  but  results  in  a  distribution  of  noise 
about  the  zero  line.  Also  plotted  in  this  figure  is  the  convolution  of  the  E- 
pulse  with  a  noisy  waveform  representing  a  target  5%  longer.  It  is  quite  easy 
to  separate  the  effects  of  noise  and  target  length  sensitivity,  suggesting 
that  random  noise  will  not  interfere  with  target  discrimination. 
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VI.  Experimental  Verification  of  the  E-pulse  Concept 

Time  domain  measurements  of  complex  conducting  target  responses  provide 
the  means  for  a  practical  test  of  the  E-pulse  concept.  The  present  experiment 
involves  measuring  the  near  scattered  field  response  of  a  simplified  aircraft 
model  to  transient  pulse  excitation.  A  Tektronix  109  pulse  generator  is  used 
to  provide  a  quasi-rectangul  ar  400  V  incident  pulse  of  nanosecond  duration. 
Transmission  of  the  pulse  over  a  5  x  6  m  conducting  ground  screen  is 
accomplished  by  using  an  (imaged)  biconical  antenna  of  axial  height  2.5  m, 
half  angle  of  8°  and  characteristic  impedance  of  160q, while  reception  is 
implemented  using  a  short  monopole  E-field  probe  of  length  1.6  cm.  Although 
the  receiving  probe  is  not  positioned  in  the  far  field  region  of  the 
scatterer,  the  resulting  measurements  have  the  desired  modal  content  in  the 
late-time  period,  provided  the  small  probe  purel y  different!’ ates  the  waveform. 
Lastly,  discrete  sampling  of  the  time  domain  waveform  is  accomplished  by  using 
a  Tektronix  sampling  oscilloscope  (S2  sampling  heads,  75  psec  risetime) 
coupled  to  a  Radio  Shack  model  III  microcomputer . 

In  this  experiment,  an  attempt  is  made  to  discriminate  between  two 
aircraft  models  by  employing  the  E-pulse  technique.  Figures  11  and  12  show 
the  measured  pulse  responses  of  simplified  Boeing  707  and  McDonnel  Douglas  F- 
18  aircraft  models,  respectively.  Each  model  is  constructed  of  aluminum  and 
has  a  geometry  as  indicated  in  the  figures.  Also  shown  in  the  figures  are  the 
dominant  natural  frequencies  extracted  from  the  late-time  portion  of  the 
response  using  a  nonlinear  least-squares  curve  fitting  technique  [13].  E- 
pulse  waveforms  can  then  be  constructed  to  atinbl  these  frequencies. 

Pulse  function  based  natural  E-pulses  constructed  to  annul  each  of  the 
two  target  responses  are  shown  in  Figures  13  and  14.  Discrimination  between 
the  targets  is  accomplished  by  convolving  these  waveforms  with  the  measured 
responses  of  Figures  11  and  12.  Figure  15  shows  the  convolution  of  the  F-18  E- 
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pulse  with  the  F  -  IS  measured  response.  Compared  to  early-time,  the  late-time 
region  has  been  effectively  annulled.  In  contrast,  Figure  16  shows  the 
convolution  of  the  707  E-pulse  with  the  F-18  measured  response.  The  result  is 
a  relatively  larger  late-time  amplitude.  Similarly,  Figure  17  displays  the 
convolution  of  the  707  E-pulse  with  the  707  measured  response.  Again,  the 
late-time  region  of  the  convolution  exhibits  small  amplitude.  Lastly,  Figure 
18  shows  the  convolution  of  the  F-18  E-pulse  and  the  measured  response  of  the 
707  model.  As  before,  the  "wrong"  target  is  exposed  by  its  larger  late-time 
convolution  response. 

VII.  Summary  and  Conclusion 

Radar  target  discrimination  based  on  the  natural  frequencies  of  a 
conducting  target  has  been  investigated.  The  response  of  targets  to  a 
particular  class  of  waveforms  known  as  "E-pulses"  has  been  demonstrated  to 
provide  an  effective  means  for  implementing  a  discrimination  process  in  the 
presence  of  random  noise. 

Two  types  of  E-pulses  have  been  identified,  natural  and  forced. 
Discrimination  based  on  natural  E-pulses  and  the  response  of  a  thin  cylinder 
target  has  been  demonstrated  theoretically. 

Most  encouraging  are  the  experimental  results  which  reveal  that  two  quite 
complicated,  similar  sized  targets  can  be  adequately  and  convincingly 
discriminated  using  natural  E-pul se  waveforms.  Further  experimentation  using 
more  accurate  aircraft  models  is  currently  being  undertaken. 
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L/a  =  200 


1  -0.0828  +  j  0 . 9  2  5 1 

2  -0.121:  jl  -012 

3  -0. 1491  +  j2.884 
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9  -0.2522  +  J8.807 

10  -0.2648  +  J9.800 


Figure  3  Orientation  for  thin  cylinder  excitation  and  first  ten  natural 
frequencies. 


>ul$e  function  based  £-pulse  synthesized 
lodes  of  thin  cylinder  target. 

>)  Natural  E-pulse  compared  to  Kennaugh 
b)  Forced  E-pulse  of  duration  Te=2.3  L/ 


Figure  7  Natural  E-pulse  constructed  using  (a)  Fourier  cosine  basis  functions 
to  eliminate  5,  6,  7,  8,  9,  and  10  modes;  (b)  pulse  function  basis 
functions  to  eliminate  3,  5,  7,  and  10  modes. 
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Figure  8  Convolution  of  10  mode  natural  thin  cylinder  E-pulse  with  60°  thin 
cylinder  impulse  response  and  60°  response  of  a  cylinder  5%  longer. 
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Figure  9  Thin  cylinder  60°  impulse  response  generated  from  first  ten  natural 
modes,  with  10%  random  noise  added. 
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Figure  11  Measured  response  of  a  Boeing  707  aircraft  model  and  seven  dominant 


natural  frequencies. 
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Figure  12  Measured  response  of  a  McOonnel  Douglas  F-18  aircraft  model  and  five 
dominant  natural  frequencies. 
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2.3  FREQUENCY  DOMAIN  E-PULSE  SYNTHESIS  AND  TARGET  DISCRIMINATION 


ABSTRACT 


A  frequency  domain  approach  to  the  B-pulse  radar  target  discrimination 
scheme  is  introduced.  This  approach  is  shown  to  allow  easier  interpretation 
of  E- pulse  convolutions  via  the  E-pulse  spectrum,  and  leads  to  a  simplified 
calculation  of  pulse  basis  function  amplitudes  in  the  E-pulse  expansion. 
Experimental  evidence  obtained  using  aircraft  models  verifies  the  single-mode 
discrimination  scheme,  as  well  as  the  aspect- independent  nature  of  the  E-pulse 
technique.  This  leads  to  an  integrated  technique  for  target  discrimination 
combining  the  E-pulse  with  single  mode  extraction  waveforms. 
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I.  INTRODUCTION 


The  time-domain  scattered  field  response  of  a  conducting  target  has  been 
observed  to  be  composed  of  an  early-time  forced  period,  when  the  excitation 
field  is  interacting  with  the  scatterer,  followed  immediately  by  a  late-time 
period  during  which  the  target  oscillates  freely  [1],  [2].  The  late-time 
portion  can  be  decomposed  into  a  finite  sum  of  damped  sinusoids  (excited  by  an 
incident  field  waveform  of  finite  usable  bandwidth),  oscillating  at 
frequencies  determined  purely  by  target  geometry.  The  natural  resonance 
behavior  of  the  late-time  portion  of  the  scattered  field  response  can  be 
utilized  to  provide  an  aspect- independent  means  for  radar  target  discrimina¬ 
tion  [3-7]. 

An  extinction  pulse  (E-pulse)  is  defined  as  a  finite  duration  waveform 
which,  upon  interaction  with  a  particular  target,  eliminates  a  pre-selected 
portion  of  the  target's  natural  mode  spectrum.  By  basing  E-pulse  synthesis  on 
the  target  natural  frequencies,  the  E-pulse  waveform  is  made  aspect- 
independent. 

Discrimination  is  accomplished  by  convolving  an  E-pulse  waveform  with  the 
measured  late-time  scattered  field  response  of  a  target.  If  the  scattered 
field  is  from  the  anticipated  target,  the  convolved  response  will  be  composed 
of  an  easily  interpreted  portion  of  the  expected  natural  mode  spectrum.  If 
the  target  is  different  from  that  expected,  a  portion  of  its  dissimilar 
natural  mode  spectrum  will  be  extracted,  resulting  in  an  unexpected  convolved 
response. 

This  paper  will  consider  the  construction  of  two  types  of  E-pulse  wave¬ 
forms.  The  first  is  designed  to  eliminate  the  entire  finite  expected  natural 
mode  spectrum  of  the  target,  and  the  second  to  extract  just  a  single  mode. 
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Both  time  and  frequency  domain  analyses  will  be  used  to  develop  a  set  of 
defining  E-pulse  equations.  The  frequency  domain  approach  will  be  shown  to  be 
of  significant  value  not  only  for  providing  a  much  more  convenient  way  of 
viewing  E-pulse  discrimination,  but  also  for  constructing  certain  E-pulse 
waveforms. 

Lastly,  experimental  results  using  aircraft  models  will  be  presented, 
demonstrating  single  mode  discrimination  and  verifying  the  aspect  independence 
of  the  E- pulse  technique. 


II.  TIME  DOMAIN  SYNTHESIS 


Assume  that  the  measured  time  domain  scattered  field  response  waveform  of 
a  conducting  radar  target  can  be  written  during  the  late-time  period  as  a 
finite  sum  of  damped  sinusoids 


r  (t)  = 


N 

V" 

n=l 


ant 


an  e 


cos  ( -nt  + 


V 


t  >  Tc 


(1) 


where  an  and  are  the  aspect  dependent  amplitude  and  phase  of  the  n'th 
mode,  T  ^  describes  the  beginning  of  late-time,  and  only  N  modes  are  assumed 
to  be  excited  by  the  incident  field  waveform.  Then,  the  convolution  of  an  E- 
pulse  waveform  e(t)  with  the  measured  response  waveform  becomes 


c(t)  =  e(t)*r(t)  = 


e(t')r  (t-t')dt' 


*  Z  ane  [V08^  +  V  +  +  >  tL 

n=l  u  J 


Tl  +  Te 


(2) 


58 


where 


and  Te  is  the  finite  duration  of  e(t) . 

Two  interesting  waveforms  are  now  considered.  Constructing  e(t)  to 
result  in  c(t)=0,  t>TL,  requires 


^  ^  =  0  l<nlN  (4) 

This  approach  was  considered  in  [3].  In  addition,  e(t)  can  also  be 
constructed  so  that  c(t)  is  composed  of  just  a  single  mode.  In  this  case  e(t) 
is  termed  a  "single  mode  extraction  signal,"  as  discussed  in  1 5—7 ] .  If  the 
specific  value  of  the  phase  constant  of  c(t)  is  unimportant,  e(t)  can  be 
synthesized  by  demanding 


Ah  =  Bn  =  0  l<nlN,  rt^m 

to  excite  the  m'  th  natural  mode.  On  the  other  hand,  requiring 


(5) 


sn=Bn  =  ° 
*m  “  0 


l<n<N,  n^m 


results  in 


am^ 


c(t)  =  cs(t)  «  ^  sin  (umt  +  4>m) 


(6) 


(7) 


and  requiring 


yields 


Ah  =  Bh  =  0  l£n£N,  n^m 


Bm  =  0 


c(t)  =  cc(t)  =  ajjj  e  cos(^mt  +  *m) 


(8) 


(9) 


The  E-pulse  resulting  from  (5)  is  termed  a  "sin/ cos"  single  mode  extraction 
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waveform/  since  c(t)  contains  both  sine  and  cosine  components.  Similarly/  (6) 
results  in  a  "sine"  and  (8)  a  "cosine"  single  mode  extraction  waveform.  With 
the  proper  normalization  of  e(t)  (giving  -  Bm),  the  convolved  waveforms  (7) 
and  (9)  can  be  combined  to  yield  plots  of  the  m' th  mode  frequencies  versus 
time,  via 


m 


t  + 


m 


tan 


-1 


Vt) 

cc(t) 


~mt  +  lo<3  IVtal  =Tl0<3 


(10) 


III.  FREQUENCY  DOMAIN  SYNTHESIS 

The  convolution  of  the  E- pulse  waveform  with  the  measured  response  wave¬ 
form  can  also  be  written  in  the  form 


N 


-nt 


c(t)  ®n  |EM  e  '  cos(V+  V  t>TI 

n=l 


where 


E(s) 


e(t)  = 


L 


i  -st 
e(t)e  dt 


is  the  Laplace  transform  of  the  B-pulse  waveform,  and 


n  =  *n  +  I- 


E, 


m 


rn 


where 


Ein  =  Im 


E(sn) 


V  / 


Ern  =  E^ 


(11) 


(12) 


(13) 


(14) 


6o 


l<n£N 


(15) 


Now»  c(t)=0  for  t>TL  requires 


Ein 


Jrn 


=  0 


or,  equivalently, 

E(sn)  =  E(sn*)  =  0  l<n^N  (16) 

In  addition,  a  sin/cos  m'th  mode  extraction  waveform  can  be  synthesized  via 

=  E(sn*)  =  0  Un^N,  n*n  (17) 

while  a  sine  m' th  mode  extraction  waveform  requires 

E(sn)  =  E(sn*)  =  0  l<n^N,  r^m  (18) 

=-E'V) 

and  a  cosine  m'th  mode  extraction  waveform  necessitates 


E(sn)  =  E(sn*)  =  0  l$n<N,  r#m 


(19) 


It  is  easily  shown  that  the  frequency  domain  and  time  domain  requirements 
for  synthesizing  an  E-pulse  are  identical.  By  expanding  the  exponential  in 
(12),  one  can  show  that  (17),  (18),  and  (19)  are  equivalent  to  (5),  (6),  and 
(8),  respectively. 

A  significant  benefit  of  using  a  frequency  domain  approach  comes  via  the 
increased  intuition  allowed  by  equation  (11).  When  an  E-pulse  waveform  is 
convolved  with  the  measured  response  of  an  unexpected  target,  the  amplitudes 
of  the  resulting  natural  mode  components  are  determined  by  evaluating  the 
magnitude  of  the  spectrum  of  e(t)  at  the  natural  frequencies  of  the  target  (a 
result  of  the  Cauchy  residue  theorem).  Thus,  the  E-pulse  spectrum  becomes 
the  key  tool  in  predicting  the  success  of  E  pulse  discrimination. 


IV.  E-PULSE  SYNTHESIS  USING  FREQUENCY  DOMAIN  APPROACH 


To  implement  the  E-pulse  requirements,  it  becomes  necessary  to  represent 
the  waveform  mathematically.  Let  e(t)  be  composed  of  two  components 

e (t)  =  ef(t)  +  ee(t)  (20) 

Here  ef(t)  is  a  forcing  component  which  excites  the  target,  and  ee(t)  is  an 
extinction  component  which  extinguishes  the  response  due  to  ef(t).  The 
forcing  component  is  chosen  freely,  while  the  extinction  component  is  deter¬ 
mined  by  first  expanding  in  a  set  of  basis  functions 

M 

ee<t>  *  E  Vm<«  (21) 

m=l 

and  then  employing  the  appropriate  E-pulse  conditions  from  section  III.  For 
an  E-pulse  designed  to  extinguish  all  of  the  modes  of  the  measured  response, 
using  (16)  results  in  the  matrix  equation 
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where 


Fm  <»>  *  of  |«m  tW} 

Ef  (s)  =o^fef  (t)J 


(23) 


and  M=2N  is  chosen  to  make  the  matrix  square.  Similar  equations  can  be 
constructed  to  accomodate  the  requirements  given  by  (17),  (18)  or  (19). 

As  in  [3] ,  two  types  of  E-pulses  can  be  easily  identified.  When  ef(t)/0, 
the  forcing  vector  on  the  right  hand  side  of  (22)  is  nonzero,  and  solutions 
for  the  basis  function  amplitudes  exist  for  any  choice  of  E-pulse  duration, 
Te,  which  does  not  cause  the  matrix  to  be  singular.  In  contrast,  when  ef(t)=0 
the  matrix  equation  becomes  homogeneous,  and  solutions  for  ee(t)  exist  only 
for  specific  durations  Te,  which  are  calculated  by  solving  for  the  zeroes  of 
the  determinantal  equation.  The  former  type  of  E-pulse  is  termed  "forced"  and 
the  latter  "natural."  Since  a  natural  E-pulse  has  no  forcing  component,  it  is 
viewed  as  extinguishing  its  own  excited  response. 

The  frequency  domain  approach  makes  it  possible  to  visualize  an  improved 
E-pulse  waveform  whose  spectrum  has  been  shaped  to  accentuate  the  response  of 
a  known  target.  For  example,  by  using  damped  sinusoids  or  Fourier  cosines  as 
basis  functions  in  the  E-pulse  expansion,  it  is  possible  to  concentrate  the 
energy  of  the  E-pulse  near  pre— chosen  frequencies,  and  enhance  the  single  mode 
response  of  a  particular  target  [4]. 


V.  CALCULATION  OF  PLLSE  BASIS  FUNCTION  AMPLITUDES 

A  very  useful  application  of  the  frequency  domain  approach  results  from 
using  subsectional  basis  functions  in  (21).  Let 


(nt-1)  A  <  t  i  m  A 
elsewhere 


(24) 
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where  g(t)  is  an  arbitrary  (but  Laplace  transformable)  function,  and  I  is  the 


pulse  width.  Then 


Fm  - 


g(t-[m-lj A )  e-st  dt 


(25) 


=  F^s)  esAe-smA 


and  the  matrix  equation  (22)  can  be  written  for  the  case  of  the  natural  E- 


pulse  as  | — 
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(27) 


Equation  (26)  is  homogeneous,  and  thus  has  solutions  only  when  the 
determinant  of  the  matrix  is  zero.  As  the  determinant  is  of  the  Vandermonde 
type  [8],  the  condition  for  a  singular  matrix  can  be  calculated  easily  as 


pir 

A  =  -  p*l,2,3. . .  l£k£N  (28) 

uk 


Thus,  the  duration  of  the  E-pulse  depends  only  on  the  imaginary  part  of  one  of 
the  natural  frequencies.  With  A  determined,  the  basis  function  amplitudes 
can  be  calculated  using  Cramer's  rule  and  the  theory  of  determinants  [9]  as 
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m  ■  l-l!"1  P(2N-1)  -  (m-i)  <29» 

where  P  is  the  sum  of  the  products  n-i  at  a  time,  without  repetitions,  of 
the  quantities  z±,  z±*,  z^^-t  zN*. 

Note  that  g(t)  does  not  appear  in  this  analysis,  and  thus  the  resulting 
pulse  amplitudes  are  independent  of  the  individual  pulse  shapes.  However, 
when  discriminating  between  different  targets,  g(t)  manifests  itself  through 
the  term  F-^s). 

VI.  EXPERIMENTAL  RESLLTS 

This  section  will  address  two  important  topics.  First,  discrimination 
between  two  similar  aircraft  models  will  be  demonstrated  using  single  mode 
extraction  waveforms.  Second,  the  aspect  independence  of  the  E-pulse 
technique  will  be  confirmed  using  the  same  two  models. 

The  measured  near  field  responses  of  simplified  McDonnel  Douglas  F-18 
and  Boeing  707  aircraft  models  have  been  published  previously  [3J.  The  models 
are  constructed  of  aluminum,  and  are  shown  in  Fig.  1.  The  dominant  natural 
frequencies  have  been  extracted  from  the  late- time  portions  of  the  responses 
using  a  continuation  method  [10].  With  these,  pulse  function  based  natural 
sine  and  cosine  single  mode  extraction  waveforms  of  minimum  duration  can  be 
constructed  using  equations  (18)  and  (19)  together  with  (24)  to  extract  the 
first  and  fourth  modes  of  the  707.  The  first  mode  waveforms  are  shown  in 
Figure  2.  Discrimination  between  the  F-18  and  the  707  can  be  accomplished  by 
convolving  the  E-pulse  waveforms  with  each  of  the  measured  aircraft  responses. 
If  the  E-pulses  are  convolved  with  the  expected  (707)  response,  the  result 
should  be  either  a  pure  first  or  fourth  mode  damped  sinusoid.  If  the  E-pulses 
are  convolved  with  the  unexpected  (F-18)  response,  the  result  will  be  an 
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unrecognizable  conglomeration  of  the  modes  of  the  unexpected  target. 

Figures  3-6  show  the  results  of  convolving  the  707  E-pulses  with  both  the 
707  and  F-18  responses.  These  frequency  plots  are  obtained  from  the  actual 
convolved  waveforms  by  using  equation  (10).  The  dotted  lines  represent  the 
slopes  of  the  expected  first  or  fourth  mode  frequencies.  It  is  seen  that  the 
frequency  plots  for  the  expected  target  (Figs  3  and  4)  parallel  the  expected 
frequency  lines  in  the  late-time,  while  those  for  the  unexpected  target  (Figs 
5  and  6)  do  not.  Thus,  the  707  and  the  F-18  are  easily  discriminated. 

It  is  also  quite  important  to  verify  experimentally  the  aspect  indepen¬ 
dence  of  the  E-pulse  technique.  To  accomplish  this,  measurements  have  been 
made  of  the  near  field  pulse  response  of  the  707  and  F-18  models  with  the 
fuselage  axes  aligned  at  various  angles  with  respect  to  the  transmitting/ 
receiving  antenna  configuration.  Two  natural  rectangular  pulse  based  E-pulses 
of  minimum  duration,  one  designed  to  eliminate  all  of  the  modes  detected  in 
the  707  responses  and  one  designed  to  eliminate  all  the  modes  in  the  F-18 
responses,  are  then  constructed  according  to  equations  (28)  and  (29). 

Convolution  of  the  707  and  F-18  E-pulses  with  each  of  the  measured  air¬ 
craft  responses  yields  the  results  shown  in  Figs  7-10.  It  is  apparent  from 
Fig.  7  that  convolving  the  F-18  E-pulse  with  the  measured  F-18  response  gives 
a  waveform  with  negligible  amplitude  in  the  late-time,  regardless  of  the 
aspect  angle.  Similar  results  are  obtained  in  Figure  8  when  the  707  E-pulse 
is  convolved  with  the  707  measured  response.  In  contrast,  Figs.  9  and  10 
reveal  that  when  the  707  E-pulse  is  convolved  with  any  F-18  response,  or  when 
the  F-18  E-pulse  is  convolved  with  any  707  response,  the  late-time  portion  of 
the  resulting  waveform  has  significant  amplitude.  Thus,  discrimination 
between  the  two  aircraft  models  is  possible  regardless  of  the  target  aspect. 

The  experimental  results  suggest  the  feasibility  of  an  E-pulse  target 


discrimination  scheme  which  utilizes  both  a  waveform  designed  to  eliminate  all 
the  modes  of  a  target  and  a  set  of  waveforms  designed  to  extract  various 
individual  target  modes.  This  technique  integrates  the  single  mode  extraction 
concept  and  the  usual  E-pulse  technique  into  a  scheme  which  has  a  greater 
potential  for  accurate  target  discrimination  decisions. 

VII.  CONCLUSIONS 

The  E-pulse  radar  target  discrimination  concept  has  been  expanded  upon, 
incorporating  single  mode  extraction  waveforms  into  an  integrated  technique. 
Both  time  and  frequency  domain  analyses  have  been  included,  and  their  equiva¬ 
lence  demonstrated.  The  frequency  domain  viewpoint  is  helpful  for 
interpreting  E-pulse  discrimination,  and  it  has  been  applied  to  pulse  basis 
function  amplitude  calculation. 

Experimental  results  obtained  using  aircraft  models  have  demonstrated  the 
validity  of  radar  target  discrimination  based  on  single  mode  extraction  wave¬ 
forms.  Most  importantly,  experimental  evidence  has  also  shown  the  E-pulse 
technique  to  be  successful  regardless  of  target  aspect  angle.  This  aspect- 
independence  of  the  E-pulse  concept  is  fundamental  for  its  application  to 
practical  situations. 
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Fig.  1.  Aluminum  aircraft  models  of  Boeing  707  (left)  and  McDonnel- 
Douglas  F-18  used  in  the  experiment.  Note  that  the  models 
are  not  to  the  same  scale,  but  rather  are  of  similar  physical 
size. 
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Fig.  2.  Natural  rectangular  pulse  function  based  first  mode  extraction 

waveforms  for  the  Boeing  707  aircraft  model. 
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response  of 


Single  mode  angular  frequency  and  damping  coefficient  extracted  from  the  convolved  outputs  of 
the  fourth  mode  extraction  signals  for  the  707  target  and  the  707  measured  response. 
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the  first  mode  extraction  signals  for  the  707  target  and  the  F-18  measured  response. 


response  of 


relative  amplitude  relative  amplitude 


■ -pul se]*[F-18  pulse  response  ( »  ■  0  j  ]  °!  [f-13  t-oul se]*[e-l 9  Pulse  response  (»  *  303)] 


time  (ns)  time  (ns) 


tim*  (ns)  time  (ns) 

Fig.  7  Convolved  outputs  of  F-18  E-pulse  and  pulse  refponses  of  F-18  model 

measured  at  (a)  0°  aspect  angle,  (b)  30°  aspect  angle,  (c)  45°  aspect 
angle,  (d)  60°  aspect  angle. 
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2.4  EXTRACTION  OF  THE  NATURAL  FREQUENCIES  OF  A  RADAR  TARGET  ROM  A  MEASURED 


RESPONSE  USING  E-PULSE  TECHNIQUE 


ABSTRACT 

A  new  scheme  is  introduced  for  extracting  the  natural  resonance 
frequencies  of  a  radar  target  from  a  measured  response.  The  method  is  based 
on  the  E-pulse  technique  and  is  shown  to  be  relatively  insensitive  to  random 
noise  and  to  estimates  of  modal  content.  Verification  of  the  technique  is 
accomplished  by  comparing  the  natural  frequencies  extracted  from  the  measured 
responses  of  a  thin  cylinder  and  a  circular  loop  with  those  obtained  from 
theory.  The  applicability  of  the  technique  to  low-Q  targets  is  also 
demonstrated,  using  the  measured  response  of  a  scale  model  aircraft. 
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I.  INTRODUCTION 


Many  recent  radar  target  discrimination  schemes  have  utilized  the  late¬ 
time  natural  oscillation  beha/ior  of  conducting  targets  [1-5].  These 
techniques  are  based  upon  che  assumption  that  the  late- time  scattered  field 
response  of  the  target  obeys  the  natural  mode  representation 

N  ~  t 

Es(t)  =  5Z  ane  n  cos ( +  -n*  t  >  TL  (1) 

n=l 

where  sn  =  an  +  j^n  is  the  aspect  independent  natural  frequency  of  the  n'th 
target  mode,  an  and  i>  n  are  the  aspect  and  excitation  dependent  amplitude  and 
phase  of  the  n'th  target  mode,  T^  describes  the  beginning  of  the  late-time 
period,  and  the  number  of  modes  in  the  response,  N,  is  determined  by  the 
finite  frequency  content  of  the  waveform  exciting  the  target. 

T.nce  the  natural  frequencies  of  the  target  are  aspect  independent,  they 
form  an  ideal  set  of  discriminants.  Employing  this  set  requires  the 
knowledge  of  the  natural  frequencies  of  a  wide  variety  of  targets.  As  a 
theoretical  determination  of  the  natural  frequencies  of  a  complex  target  is 
impractical,  it  becomes  necessary  to  develop  a  scneme  for  extracting  the 
frequencies  from  a  measurement  of  the  response  of  the  target. 

A  typical  approach  might  use  Prony's  method  [6,7].  Although  this 
technique  is  simple  and  efficient,  it  has  been  found  to  be  overly  sensitive 
to  random  noise  and  to  the  number  of  modes  assumed  to  be  present  in  the 
measured  response.  A  nonlinear  ]east  squares  curve-fitting  scheme  can 
overcome  these  drawbacks,  but  requires  a  time  consuming  minimization  involving 
4N  variables  [8] . 

In  this  paper  an  alternative  scheme  is  proposed,  utilizing  E-pulse 


waveforms.  It  will  be  shown  that  this  approach  also  overcomes  the  drawbacks 
of  Prony's  method,  but  is  more  efficient  than  nonlinear  curve-fitting.  The 
scheme  was  first  introduced  in  [12]  and  [15],  and  is  based  on  the  frequency 
domain  E-pulse  concept  first  discussed  in  [14]. 

II.  THE  E-PtLSE  TECHNIQUE 

An  extinction-pulse  (E-pulse)  e(t)  is  defined  as  a  waveform  of  finite 
duration  T0  which  extinguishes  Es(t)  in  the  late-time  [1].  That  is, 
convolution  of  e(t)  and  £s(t)  yields  the  null  result 

Je 

c(t)  =  e(t)*Es(t)  =  I  e(t')Es(t-t')dt'  =  0  t  >  TL  +  T0  (2) 

o 

□sing  the  natural  mode  representation  for  Es(t)  allows  (2)  to  be  written  as 


c(t)  =  Z  l5^)}  ane  n  cos(uJnt  + 
n=l 

where  in  is  a  new  phase  factor  and  E(s) 


V  -  0  t  >  TL  +  Te  (3) 

is  the  Laplace  transform  of  e ( t ) 


E(s) 


(4) 


If  the  natural  frequencies  of  a  target  are  known,  an  E-pulse  for  that  target 
can  be  synthesized  by  denanding 

E(sn)  =  E(sn*)  =  0  1  ^  n  S  N  (5) 

and  the  convolution  (3)  will  yield  zero  regardless  of  the  aspect  angle  of  the 
target  for  which  Es(t)  is  measured. 


Conversely,  if  the  natural  frequencies  of  a  target  are  unknown,  they  can 
be  extracted  from  the  measured  Es(t)  by  solving  the  integral  equation  given  by 
(2).  This  equation  can  either  be  solved  directly  for  the  complex  frequencies 


used  to  construct  e(t)  or  for  e(t)  itself.  If  e(t)  is  determined,  the  complex 
frequencies  eliminated  by  e(t)  can  be  found  by  locating  the  roots  to  [12] 


E(s)  =  0 


(6) 


An  E-pulse  waveform  can  be  represented  as  an  expansion  over  a  set  of 
basis  functions  jf^Ct)  j 


K 

e(t)  =  £skfk  (t)  (7) 

fel  K  K 


where  ak  are  the  basis  function  amplitudes.  Then,  the  requirements  (5) 
result  in  a  matrix  equation  for  the  real  amplitude  parameters  a^.  Choosing 
K=2N  results  in  a  homogeneous  equation  which  tes  solutions  only  at  discrete 
values  of  Tg. 

The  difficulty  associated  with  locating  all  the  roots  of  (6)  can  be 
overcome  by  choosing  subsectional  basis  functions  in  the  E-pulse  expansion 


f„(t)  = 


(8) 


g(t  -  [k-lj  A  )  (k-l)A  <  t  <  k  A 

0  elsewhere 

resulting  in  the  E-pulse  shown  in  figure  1.  Here  g(t)  is  any  Laplace- 
transformable  function.  The  E-pulse  spectrum  is  then  easily  calculated  as 


S  A 


E  (s)  =  Fj_(s)  e  S  a  i,e 


-sk  A 


(9) 


k=l 


where  F1(s)  is  the  Laplace  transform  of  the  first  basis  function.  Now  the 
roots  to  E (s)  are  easily  found  by  solving  the  polynomial  equation 


K 

Lt  =  0  (io> 

k=l  K 

where 

-SA 

Z  =  e  (11) 


This  set  of  functions  also  allows  a  simple  calculation  of  duration  Tg  when 
synthesizing  E-pulses.  It  has  been  shown  [9]  that  the  solutions  to  the 
homogeneous  matrix  equation  demanded  by  (5)  are 

i=PJL  1  ^  k  <  N,  p  =  1,2,3,...  (12) 

~k 


III.  SOLUTIONS  TO  THE  INTEGRAL  EQUATION 

Solutions  to  the  integral  equation  (2)  have  been  obtained  using  two 
different  methods.  The  first  approach  is  to  minimize 

c2 ( t )  =  [e (t)  *  Es(t)J2  (13) 

over  the  range  of  t  corresponding  to  the  late-time  period  of  the  convolution, 
t>T^+Te.  Hiis  can  be  done  either  with  respect  to  the  complex  frequencies  used 
to  construct  the  E-pulse  via  (5),  or  with  respect  to  the  basis  function 
amplitudes  used  to  construct  e(t)  via  (7).  If  minimization  is  done  with 
respect  to  the  complex  frequencies,  then  the  estimated  values  of  the  natural 
frequencies  contained  in  Es(t)  are  available  directly  at  the  minimum  point. 
If  minimization  is  done  with  respect  to  the  basis  function  amplitudes,  then 
estimates  of  the  natural  frequencies  in  Es(t)  are  obtained  by  solving  (6)  or 
(10)  using  the  values  of  obtained  at  the  minimum  point.  The  latter 
minimization  is  computationally  less  intensive,  since  a  new  E-pulse  need  not 
be  created  at  each  step  in  the  minimization,  but  it  requires  a  normalization 
scheme  to  prevent  the  trivial  solution  a  ^ 

The  integral  equation  can  also  be  bv^.j  by  using  the  method  of  moments 
[10].  The  E-pulse  is  expanded  as  in  (7)  and  moments  are  taken  with  a  set  of 
weighcing  functions  |wm(t)\ 


K  Te 

<^m(t)'  £  \  /  fk(t')Es(t-t,)dt,)>  =  0  m=l,2,...,M 

\  k=l  o 
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where  the  brackets  <  >  indicate  the  inner  product 

<f(t),  g(t)>  =  /f(t)g(t)dt  (15) 

Choosing  K=M=2N  to  reflect  the  number  of  modes  believed  to  be  in  Es(t) 
results  in  a  homogeneous  matrix  equation  for  the  basis  function  amplitudes,  a 
solution  to  which  is  possible  only  at  discrete  values  of  Tg  which  cause  the 
determinant  of  the  matrix  of  coefficients  to  vanish.  By  choosing  subsectional 
basis  functions  (8)  an  estimate  for  the  natural  frequencies  in  Es(t)  is  found 
by  solving  (10). 

With  either  of  the  two  approaches,  the  amplitudes  and  phases  of  the 
natural  modes  comprising  Es(t)  can  be  found  as  a  last  step  by  using  linear 
least  squares. 

IV.  DISCUSSION  OF  THE  METHODS 

The  major  benefit  of  each  of  the  two  E-pulse  methods  is  that  they 
overcome  the  sensitivity  of  Prony's  method  to  noise  and  to  estimates  of  modal 
content.  Although  the  minimization  approach  requires  a  computationally  time 
consuming  multi-variable  minimization  procedure,  because  the  amplitudes  and 
phases  are  not  involved  the  total  number  of  parameters  utilized  is  only  2N, 
compared  to  4N  for  nonlinear  curve-filling,  and  only  half  the  number  of 
initial  guesses  are  needed  for  starting  the  minimization.  On  the  other  hand, 
the  moment  method  approach  is  extremely  efficient;  no  minimization  scheme  is 
required,  only  a  search  for  the  single  parameter  Tg. 

The  insensitivity  of  these  two  approaches  to  the  number  of  modes 
estimated  to  be  present  in  Es(t)  leads  to  a  particularly  useful  scheme. 
Extraction  of  the  frequencies  in  Es(t)  may  begin  by  assuming  a  small  number  of 
modes  (usually  one)  to  be  present.  For  the  methods  requiring  minimization, 
the  frequencies  extracted  with  a  small  number  of  modes  assumed  present  prove 


to  be  very  good  initial  guesses  for  the  case  when  more  modes  are  assumed 
present.  For  the  moment  method,  because  of  the  relationship  (12),  the  value 
of  Te  obtained  allows  an  excellent  estimate  of  the  Tg  to  expect  with  more 
modes  allowed. 

The  reasons  for  the  success  of  the  E-pulse  method  in  the  presence  of 
random  noise  can  be  identified  by  investigating  the  relationship  between 
Prony's  method  and  the  moment  method  approach  to  the  E-pulse  technique.  By 
choosing  impulses  for  both  expansion  and  weighting  functions  in  (14)  and 
taking  K=M+1,  an  inhomogeneous  matrix  equation  results.  The  matrix  can  be 
solved  directly  for  ak  and  the  natural  frequencies  contained  in  Es(t) 
estimated  by  solving  (10).  Such  an  approach  is  found  to  be  identical  in  all 
respects  to  Prony's  method. 

Because  of  the  discrete  nature  of  the  convolution  integral,  each  of  the 
matrix  entries  and  the  right  hand  side  vector  elements  are  merely  a  single 
sampled  value  of  the  measured  waveform.  Thus,  Prony's  method  is  very 
sensitive  to  noise  contamination  of  Es(t)  and  so  is  an  inherently  ill- 
conditioned  algorithm  (see  [16]).  The  moment  method  approach  can  be  viewed 
as  a  generalization  of  the  basic  Prony's  method,  where  preprocessing  of  the 
measured  data  is  naturally  introduced.  The  preprocessing  step  can  be 
incorporated  in  two  places,  individually  or  simultaneously,  each  of  which 
works  to  reduce  the  noise  sensitivity  of  the  technique.  First,  by  utilizing 
expansion  functions  which  together  span  the  E-pulse  duration  Te  (such  as 
rectangular  pulses),  the  convolution  integral  in  (14)  performs  a  moving  window 
type  smoothing  of  the  measured  data.  Second,  by  using  weighting  functions 
which  together  span  the  late- time  region  of  the  convolution,  the  inner  product 
integral  (15)  introduces  additional  smoothing.  Although  the  numerical 
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examples  presented  below  utilize  impulse  functions  for  weighting,  preliminary 
results  using  rectangular  pulses  indicate  an  additional  reduction  in  noise 
sensitivity. 

Hie  result  of  preprocessing  the  measured  data  is  matrix  elements  which 
are  each  a  function  of  a  large  portion  of  the  measured  data.  Consequently, 
solutions  are  less  sensitive  to  the  perturbation  of  individual  sampled  values 
of  Es(t),  and  so  the  E-pulse  technique  is  a  better  conditioned  algorithm. 
Further  noise  reduction  capabilities  can  be  introduced  by  using  basis 
functions  that  individually  span  Te  (such  as  Fourier  cosines)  and/or  weighting 
functions  that  individually  span  the  late- time  of  the  convolution,  involving 
even  more  of  the  measured  data  in  each  matrix  element.  However,  choosing 
basis  functions  which  are  not  subsectional  complicates  the  search  for  zeroes 
of  the  E-pulse  spectrum,  since  (10)  can  no  longer  be  used. 

One  last  distinction  between  Prony’s  method  and  the  E-pulse  techniques  is 
the  manner  in  which  the  E-pulse  duration  Te  is  selected.  Whereas  in  Prongs 
method  Te  is  determined  solely  by  the  sampling  interval  (so  that  discrete 
convolution  may  be  performed),  in  the  E-pulse  technique  Te  is  linked  directly 
to  the  natural  frequencies  contained  in  the  measured  response. 

When  minimizing  with  respect  to  the  natural  frequencies,  the  duration  of 
the  E-pulse  convolved  with  the  data  is  chosen  to  be  the  smallest  allowed  by 
(12) .  Note  that  this  is  given  by 

<Te>min  -  2N  £h  <«> 

where  is  the  largest  value  of  ^  present  among  the  natural  frequencies  in 
the  response.  This  choice  results  in  the  largest  amount  of  the  late-time 
convolved  response  used  in  minimizing  c^(t).  As  the  frequencies  change 
during  the  minimization  process  the  duration  will  change  as  well,  but  it  is 
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always  related  through  (16) . 

When  minimizing  with  respect  to  the  basis  function  amplitudes,  the  E- 
pulse  duration  is  allowed  to  be  a  free  variable.  However,  at  each  step  the 
basis  function  amplitudes  represent  an  E-pulse  which  eliminates  a  certain  set 
of  natural  frequencies,  and  the  duration  is  tied  to  those  frequencies  via 
(12). 

In  employing  the  moment  method,  the  E-pulse  duration  is  the  only  free 
parameter.  When  a  zero  of  the  determinantal  equation  is  located,  it  will 
again  be  tied  to  the  natural  frequencies  eliminated  by  the  E-pulse  through 
(12).  As  before,  it  is  prudent  to  search  for  the  smallest  value  of  Te  which 
satisfies  the  determinantal  equation  since  this  results  in  the  longest  late¬ 
time  portion  of  the  convolution. 

As  a  last  note,  the  E-pulse  technique  as  a  whole  should  not  be  considered 
just  a  generalization  of  Prony's  method.  For  instance,  the  minimization  with 
respect  to  natural  frequencies  approach  completely  avoids  the  second  step  in 
Pron/s  method  —  solving  a  polynomial  equation  for  the  natural  frequencies 
(which  itself  can  be  an  ill-conditioned  problem.)  In  the  moment  method 
approach,  a  different  choice  of  basis  functions  requires  finding  not  the  roots 
of  a  polynomial,  but  those  of  a  sum  of  different  functions  altogether. 

V.  EXAMPLES 

As  a  simple  numerical  example,  Figure  2  shows  the  theoretical  impulse 
response  of  a  thin  cylinder  oriented  at  30°  with  respect  to  the  incident 
field,  constructed  using  the  first  eight  natural  frequencies  of  the  cylinder 
[11] .  An  attempt  will  be  made  to  extract  these  frequencies  from  the  response 
using  the  E- pulse  techniques. 

To  simulate  a  practical  situation,  this  waveform  is  sampled  at  500 
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points.  The  convolution  indicated  in  (2)  is  then  carried  out  by  interpolating 
linearly  between  sampled  points,  and  analytically  integrating  the  product  of 
the  linear  curve  and  the  mathematical  representation  of  the  basis  functions 
from  (8).  To  keep  the  integrals  simple,  rectangular  pulse  basis  functions  are 
used  throughout. 

The  minimization  schemes  utilize  standard  Newton's  method  and  thus 
require  initial  guesses.  For  minimization  with  respect  to  the  complex 
frequencies,  initial  guesses  for  a  n  and  w  n  are  required.  These  can  often 
be  obtained  from  the  Fourier  spectrum  of  12s  (t)  via  the  FET.  For  minimization 
with  respect  to  che  basis  function  amplitudes,  initial  guesses  are  required 
for  ak.  Anticipating  these  values  is  difficult,  since  intuition  is  lacking. 
However,  a  good  guess  for  the  amplitudes  can  be  obtained  by  constructing  at 
the  first  step  an  E-pulse  based  on  guesses  for  an  and  ^n. 

The  number  of  basis  functions  chosen  in  the  E- pulse  expansion  reflects 
the  number  of  modes  expected  —  two  functions  for  each  mode.  If  a  DC 
component  is  present  in  Es(t)  as  an  artifact  of  the  measurement  system,  it 
can  be  effectively  eliminated  by  utilizing  an  additional  basis  function  [9]. 

To  demonstrate  the  insensitivity  of  the  E- pulse  techniques  to  the  number 
of  modes  assumed  present  in  E®(t),  an  attempt  will  be  made  to  extract  four  of 
the  eight  frequencies  used  to  construct  Figure  2.  Figure  3  shows  the  results 
of  using  the  moment  method  with  impulse  weighting  functions,  and  minimizat^.n 
with  respect  to  the  complex  frequencies.  Obviously,  the  first  four  modes  have 
been  extracted  with  very  good  precision,  even  though  the  number  of  modes 
present  has  been  drastically  underestimated. 

To  demonstrate  the  insensitivity  of  the  E-pulse  techniques  to  the 
presence  of  random  noise,  an  attempt  will  be  made  to  extract  the  frequencies 
from  a  noisy  version  Figure  2.  The  amplitude  of  the  added  noise  is  chosen  to 


be  10%  of  the  maximum  value  of  the  waveform.  Figure  4  shows  the  results  of 
using  the  same  moment  method  and  minimization  approaches,  and  the  results  are 
seen  to  be  quite  adequate.  Note  that  the  values  of  ^  n  obtained  are  usually 
much  better  than  the  values  of  a  n.  As  a  direct  comparison.  Figure  3d  of  [8J 
shows  the  natural  frequencies  extracted  from  the  noisy  response  using  Prony's 
method.  It  is  apparent  that  utilizing  the  E-pulse  concept  provides  a 
decrease  in  the  noise  sensitivity  of  the  extracted  natural  frequencies. 

As  a  more  practical  example,  Figure  5  shows  the  measured  surface  current 
response  of  a  thin  cylinder  to  a  nanosecond  pulse  excitation  field  (see  [1] 
for  a  detailed  description  of  the  experiment).  The  measurement  system  has 
sampled  this  response  at  a  total  of  240  points.  Figure  6  shows  the  results  of 
extracting  the  five  dominant  natural  frequencies  using  minimization  with 
respect  to  the  complex  frequencies.  Rectangular  pulse  basis  functions  have 
been  used  again,  and  sampled  point  convolution  performed  as  described  above. 
The  comparison  of  the  extracted  frequencies  to  those  given  by  theory  is 
excellent. 

As  a  second  practical  example,  Figure  7  shows  the  measured  late-time 
scattered  field  response  of  a  thin  wire  circular  loop  to  a  nanosecond  pulse 
excitation  field.  Ihe  measurement  system  has  sampled  this  waveform  at  1024 
points.  Figure  8  shows  the  result  of  using  minimization  with  respect  to  the 
complex  frequencies  to  extract  the  six  dominant  natural  frequencies.  Also 
shown  in  Figure  8  are  the  theoretical  values  calculated  using  a  Fourier  series 
type  solution  [13].  Again,  the  frequencies  extracted  from  the  measured 
waveform  compare  quite  well  with  the  theoretical  values. 

Lastly,  it  is  important  to  demonstrate  that  the  E-pulse  technique  will 
work  successfully  for  low-Q  type  targets.  Figure  11  of  [1]  shows  the  measured 
nanosecond  pulse  scattered  field  response  of  a  Boeing  707  aircraft  model, 
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sampled  at  400  points.  Although  no  theory  is  available  for  predicting  the 
natural  frequencies  of  this  structure,  they  have  been  extracted  from  the 
measured  response  using  a  nonlinear  least  squares  technique  [8J.  Figure  9 
shows  the  six  dominant  natural  frequencies  extracted  using  the  least  squares 
approach,  and  also  those  extracted  using  minimization  of  the  late-time 
convolved  response  with  respect  to  the  complex  frequencies  (again,  using 
rectangular  pulse  expansion  functions).  The  agreement  between  the  results 
obtained  by  these  radically  different  methods  is  very  good,  verifying  the 
usefulness  of  the  E-pulse  technique  for  lcv-Q  structures. 

VI.  CONCLUSIONS 

Two  new  approaches  to  extracting  the  natural  frequencies  of  a  radar 
target  have  been  presented.  They  are  based  on  the  B-pulse  technique,  and  have 
been  shown  to  be  relatively  insensitive  to  random  noise  and  to  the  number  of 
modes  assumed  present  in  the  measured  response.  They  have  also  been  shown  to 
work  in  practical  situations,  accurately  reproducing  the  expected  natural 
frequencies  of  a  thin  cyclinder  and  a  circular  loop  from  measr  ad  values  of 
their  transient  responses.  Lastly,  the  new  methods  have  been  used  to  verify 
the  values  of  the  natural  frequencies  of  a  low-Q  aircraft  model  obtained  using 
a  different  approach. 
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Figure  1.  E-pulse  waveform  constructed  using 
subsectional  basis  functions. 
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Figure  2.  Eight  mode  30°  thin  cylinder  impulse  response. 
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Fig.  3.  Natural  frequencies  extracted  from  thin  cylinder  impulse  response, 
four  modes  assumed  present. 
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Figure  5.  Measured  surface  current  response  of  a  chin  cylinder 
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Fig.  6.  Natural  frequencies  extracted  from  surface  current  response  of 
18  1/4"  thin  cylinder. 
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Figure  8.  Natural  frequencies  extracted  from  measured 
response  of  thin  wire  circular  loop. 
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2.5  NOISE- INSENSITIVITY  OF  E- PULSE  TECHNIQUE  FOR  TARGET  C I  SC  R IM  IN  ATI  CN 


To  Drove  the  noise  insensitivity  of  our  schene,  we  have  created  very 
noisy  radar  responses  of  complex  targets  by  intentionally  adding  a  ’arne 
random  noise  to  the  measured  radar  responses  of  the  targets.  These  noisy 
responses  were  then  used  to  convolve  with  the  discriminant  signals,  the 
E-pulse  and  the  single-mode  extraction  signals,  of  the  targets.  We  have 
found  that  the  discriminant  signals  of  the  targets  are  very  powerful  and 
effective  in  rejecting  a  large  random  noise  and  are  capable  of  discriminating 
between  the  right  and  the  wrong  targets  from  very  noisy  radar  responses.  The 
following  figures  will  demonstrate  this  finding. 

Figure  1  shows  the  pulse  response  (the  response  excited  by  an  incident 
Gaussian  pulse)  of  8707  model  measured  at  90°  aspect  angle  without  adding  an 
extra  random  noise.  Figure  2  is  the  convolved  output  of  the  pulse  response  of 
Fig.  1  with  the  E-pulse  of  8707  model.  As  expected,  a  very  small  output  was 
obtained  in  the  late-time  period  of  the  convolved  output.  This  indicates  that 
the  Dulse  response  of  Fig.  1  came  from  the  right  target  of  B707  model.  Next, 
a  very  noisy  pulse  response  was  created  by  intentionally  adding  a  large 
random  noise  (created  by  a  computer)  to  the  measured  pulse  response  of  B707 
model  shown  in  Fig.  1.  This  random  noise  amounted  to  30%  of  the  maximum 
amplitude  of  the  measured  response  of  Fig.  1.  The  created  noisy  pulse 
response  is  shown  in  Fig.  3.  When  this  noisy  pulse  response  of  Fig.  3  was 
convolved  with  the  E-pulse  of  B707  model,  a  very  satisfactory  convolved 
output  was  obtained,  as  shown  in  Fig.  4.  This  convolved  output  resembles 
that  of  Fig.  2;  the  early-time  response  stayed  nearly  unchanged  and,  more 
importantly,  the  late-time  response  still  remained  small.  This  indicates  that 
the  E-pulse  of  8707  model  was  capable  of  identifying  that  the  noisy  pulse 
response  of  Fig.  4  belonged  to  8707  model. 
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Figure  1.  Pulse  response  of  B707  model  measured  at  90°  aspect 
noise  added. 
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Figure  2.  Convolved  outout  of  the  E-pulse  of  B707  model  with  the  pulse  response  of  K707 
model  measured  at  90°  aspect  without  an  extra  noise  added.  A  very  small  late¬ 
time  response  was  obtained. 
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Figure  3.  Pulse  response  of  B707  model  measured  at  90  aspect  angle  with  an  extra  random 
noise  (30%  of  the  maximum  waveform  amplitude)  added. 
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Figure  Convolved  output  of  the  E-pulse  of  B707  model  with  the  oulse  response  of  B/07 
model  measured  at  90°  aspect  angle  with  an  extra  random  noise  (30%  maximum 
amplitude)  added.  The  late-time  response  was  still  very  small. 


Next,  we  tried  to  discriminate  a  wrong  target,  a  F-13  model,  with  the 
E-pulse  of  B707  model  using  noisy  pulse  responses  of  the  wrong  target,  "igure 
5  is  the  pulse  response  or  h - 1 3  model  measured  at  90°  aspect  angle  without 
adding  a  random  noise,  when  the  response  of  Fig.  5  was  convolved  with  the 
E-pulse  of  B707,  the  convolved  output  is  shown  in  Fig.  6.  In  Fig.  6,  it  is  seen 
that  a  large  late-time  response  was  obtained.  This  indicates  that  the  resDonse 
of  Fig.  5  came  from  a  wrong  target  other  than  B707  model.  Next,  a  very  noisy 
pulse  response  of  F-18  model  was  created  by  intentionally  adding  a  larqe 
random  noise,  30%  of  the  maximum  amplitude  of  the  response  of  Fig.  5,  to  the 
measured  response  of  Fig.  5.  This  noisy  pulse  response  of  F-18  model  is 
shown  in  Fig.  7.  The  noisy  pulse  response  of  Fig.  7  was  then  used  to  convolve 
with  the  E-pulse  of  B707,  and  the  convolved  output  is  shown  in  Fig.  8.  The 
convolved  output  of  Fig.  3,  as  compared  with  the  result  of  Fig.  6,  shows  a 
relatively  unchanged  early-time  response  followed  by  a  still  large  and 
somewhat  noisy  late-time  response.  This  large  late-time  response  is  sufficient 
to  indicate  that  the  noisy  pulse  response  of  Fig.  7  belonged  to  a  wrong  target 
other  than  B707  model . 

We  have  also  convolved  noisy  pulse  responses  of  complex  targets  with 
single-mode  extraction  signals  of  the  targets.  Similar  results  as  those 
described  in  Figs.  1  to  3  were  obtained;  the  single-mode  extraction  signals 
of  a  complex  target  are  capable  of  discriminating  between  the  right  target  and 
wrong  targets  based  on  very  noisy  pulse  responses  of  the  targets. 

Based  on  this  study,  w<a  can  conclude  that  our  target  discrimination 
shceme  of  using  the  E-oulse  and  single-mode  extraction  signals  is  very 
noise  insensitive. 
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Figure  5.  Pulse  resnonse  of  F- 18  model  measured  at  90°  asoect  angle  without  an  extra 
noise  added. 
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Figure  7.  Pulse  response  of  F- 13  model  measured  at  90°  aspect 
noise  (30%  of  the  maximum  waveform  amplitude)  added 
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Figure  8.  Convolved  output  of  the  E-pulse  of  B707  model  with  the  nulse  response  of  F- 18 
model  measured  at  90^  aspect  angle  with  an  extra  random  noise  (30/6  maximum 
aolitude)  added,  The  late-time  response  remained  large. 


2.6  DETERMINATION  OF  THE  IMPULSE  RESPONSE  OF  A  COMPLEX  TARGET  THROUGH 


DECONVOLUTION 

When  a  target  is  illuminated  bv  an  electromagnetic  impulse,  the  impulse 
response  of  the  target  measured  by  a  receiving  system  is  usually  entirely 
distorted  due  to  non-ideal  frequency  characteristics  of  the  system. 
Furthermore,  the  existence  of  a  delta  function  in  the  response  due  to 
specular  reflection  from  the  target  will  exacerbate  the  problem.  We  have 
developed  a  method  which  is  capable  of  recovering  the  actual  impulse 
response  of  the  target  from  real-time  measured  data. 


(I)  Mathematical  background 

Hie  equation  we  cure  attempting  to  solve  is 

E*H  =  R  (1) 

where  E  and  R  are  measured  and  H  is  to  be  solved  for.  When  this  is 
discretized,  the  set  of  equations  for  the  values  of  E,  H  and  R  at  discrete 
time  intervals  is 

ri  *  nhi 

r2  “  elh2  +  e2hl 

r3  *  elh3  +  e2^2  +  e3hl 

rn  *  nh,  +  ••••  +  enhl 
Or,  in  matrix  form 


'  e1  0  0  .  .  0 

f  \ 

hl 

V 

e2  e!  0  .  .0 

h 2 

r2 

0 

• 

*3 

= 

r3 

;  e2  ei  o 

• 

en  .  .  ®3  e2  ei 

*n 

rn 

\  ^ 

or,  E*H  »  R  (2) 

where  E  is  now  a  matrix  and  H  and  R  are  vectors.  Here,  H  represents  the 
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impulse  response  of  a  target,  E  represents  the  system  response  of  transmission 
and  reception,  and  R  represents  measured  data  of  the  total  response.  Solving 
this  is  an  ill-conditioned  problem.  Instead,  we  minimize  the  matrix  equation 

[E-H-R]2  +  p*[H]2.  (3) 

p  is  a  parameter  which  can  be  varied  to  a  suitable  value  to  obtain  the  most 
information  with  the  least  noise,  p  =  0  corresponds  to  the  original  ill- 
conditioned  problem.  As  p  is  increased,  the  equation  becomess  less  ill- 
conditioned,  and  the  results  are  smoother.  If  p  is  made  too  large,  the 
results  are  too  smooth  and  higher  frequencies  are  lost.  Minimizing  the  matrix 
equation,  we  obtain 

Efc.E-H  +  p-H  =  Efc-R.  (4) 

This  equation  is  less  ill-conditioned  than  the  original  one.  The  product 
Efc-E,  however,  is  twice  as  ill-conditioned  as  the  original  matrix  E.  We  avoid 
forming  this  through  singular  value  decomposition.  This  is  a  method  of 
breaking  a  matrix  E  into  three  matrices  U,  Vfc  and  S  such  that 

E  =  U-S-V*1  (5) 

with  S  diagonal.  Also, 

Vfc  =  V"1 
and  Ufc=  U"1. 

Correspondingly, 

Efc  *  V* S* Ufc  (5') 

so  Efc-E  ■  V-S-U^U-S-^  *  V-S2-^ 
and  our  equation  becomes 

[V-S2-^  +  p-IJ-H  =  V-S-Ut-R.  (6) 

Vfc*V  *  I,  so 

[S2  +  p*  I]  *V^*H  -  S-U^R.  (7) 

This  equation  is  then  solved  easily  for  H  since  [S2  +  p*I]  is  diagonal  and 
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V"1  =  V*\ 


H  =  V-(S2  +  p-I)"1  S.Ufc-R.  (8) 

(II)  Practical  difficulties 

when  a  data  set  representing  a  reasonably  smooth  impulse  response  is 
convolved  with  another  data  set  representing  a  system  response,  a  third  set 
representing  measured  data  is  obtained.  When  the  system  response  is 
deconvolved  out,  the  impulse  response  is  recreated,  as  it  should  be,  with  only 
small  differences  near  the  end  of  the  data  set  due  to  data  truncation. 

If,  however,  the  impulse  response  contains  a  delta  function  within  it, 
represented  by  a  single  data  point  of  value  much  larger  than  the  surrounding 
data,  reproduction  of  the  impulse  response  is  much  more  difficult. 

The  delta  function  creates  an  instability  in  the  deconvolution  process 
separate  from  the  ill-conditioned  nature  of  the  problem.  The  source  of  the 
difficulty  lies  in  the  matrix  equation  we  are  minimizing, 

11E-H  -  R  || 2  +  P  •  1 1 H  ; | 2 .  (3) 

The  term  |]  Hj|  ,  for  the  H  which  we  are  attempting  to  discern  contains  the 
square  of  a  delta  function  as  one  of  its  terms.  When  minimization  is  carried 
out,  the  solution  will  forego  the  target  H,  (which  causes  the  first  term  in 
eq.  (3)  to  be  zero)  since  the  second  term  would  then  be  very  large.  Instead, 
another  data  set,  IT,  will  be  obtained  which  doesn't  have  a  delta  function 
term  and  for  which  the  second  term  will  be  reduced.  Obviously,  increasing  the 
parameter  p  increases  the  instability,  rather  than  lowering  it  out.  Only  as  p 
approaches  zero,  does  H'  approach  H  as  eq.  (3)  bee  ones  equivalent  to  eq.  (2). 

If  noise  is  added  to  the  data  set  representing  the  measured  response,  we 
don't  have  the  luxury  of  allowing  p  to  approach  zero.  In  order  to  obtain  a 
resonable  approximation  to  the  desired  H,  p  must  be  a  significantly  large 
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number,  or  else  the  solution  obtained  will  contain  high  frequency  oscillations 
that  overwhelm  the  desired  solution.  Thus,  attempting  to  obtain  an  impulse 
response  which  contains  a  delta  function  in  the  presence  of  noise  is 
difficult. 

This  situation  is  important  in  determining  the  impulse  response  of  a 
radar  target  since  a  delta  function  exists  in  the  response  due  to  specular 
reflection  from  the  target.  This  delta  function  response  is  there  even  though 
the  transmitted  pulse  is  band-limited. 

To  circumvent  this  difficulty,  the  delta  function  must  be  removed 
somehow.  We  can't  simply  subract  it  because  it  is  smeared  out  across  the 
measured  data.  Instead,  the  measured  data  set  is  first  convolved  with  a 
smooth  pulse  roughly  gaussian  in  shape.  This  has  the  effect  of  smearing  out 
the  delta  function  that  is  implied  in  the  measured  data  before  attempting  to 
deconvolve  the  impulse  response  out.  The  pulse  is  chosen  thin  enough  that 
the  desired  features  of  the  impulse  response  won't  be  altered  much,  but  thick 
enough  that  deconvolution  can  be  carried  out  reasonably  well. 

We  applied  this  method  to  data  simulating  an  impulse  response  containing 
a  delta  function.  We  created  a  data  set  that  represented  a  smooth  target 
impulse  response,  H,  with  a  single  point  of  large  magnitude.  We  convolved  H 
with  a  second  set  representing  the  known  system  transfer  function,  E,  to  get 
R.  We  then  added  random  noise  to  R  to  get  R*  and  attempted  to  regain  the 
target  impulse  response  by  directly  deconvolving  E  from  R'.  We  didn't 
succeed — H'  didn't  resemble  H  for  any  values  of  the  parameter  pc  We  then 
convolved  R'  with  a  thin  pulse  to  get  R",  and  then  deconvolved  E  from  this. 
Upon  doing  so,  we  obtained  a  data  set  H"  which  resembled  our  original  H, 
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except  for  some  roundedness  which  is  to  be  expected  from  convolution  with  the 
pulse. 

(Ill)  Experimental  results 

We  applied  this  method  to  actual  measured  data.  The  experimental  set-up 
consists  of  a  transmitting  antenna,  a  receiving  probe  and  a  target,  all  placed 
above  a  ground  plane.  The  direct  response  of  the  probe  to  a  pulse  without  a 
target  present  is  measured.  This  serves  as  the  reference  response,  e(t)  = 

E.  The  target  is  positioned  and  the  response  of  the  probe  to  the  scattered 
wave,  m(t)  =  M,  is  then  measured.  Since  the  probe  still  responds  to  the 
direct  wave,  the  direct  response  is  subtracted  off,  leaving  the  pure 
scattered  response,  r(t)  =  m(t)  -  e(t)  =  M  -  E  =  R.  We  then  use  (8)  via  (5) 
and  (5')  to  obtain  H  such  that,  as  near  as  possible,  E  H  =  R. 

An  example  is  given  here.  In  this  example,  the  target  was  a  6"  thin 
cylinder,  and  its  scattered  field  was  measured  by  three  different  probes  with 
the  lengths  of  2",  3"  and  4  1/2".  When  we  extracted  the  impulse  responses 
without  regularizing  the  scattered  responses,  we  obtained  three  different 
responses,  shown  in  figure  1,  with  the  4  1/2"  response  looking  especially  bad. 
we  also  extracted  the  impulse  response  after  convolving  the  scattered 
responses  with  a  smooth  pulse,  shown  in  figure  Z.  The  results  for  the  three 
cases  were  improved,  with  the  three  looking  much  more  similar  to  one  another, 
even  the  4  1/2"  case. 

3 .  Future  Plans 

The  following  topics  will  receive  major  attention  in  the  future. 

(1)  The  noise-insensitivity  of  our  scheme. 
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Fig.  1.  The  impulse  responses  of  a  6"  thin  cylinder  extracted  from  its 

scattered  fields  measured  by  three  different  receiving  probes  with 
lengths  of  2",  3"  and  4  1/2".  The  scattered  fields  were  not 
regularized  by  convolving  with  a  snooth  pulse  before  deconvolution. 
The  extracted  impulse  responses  of  the  target  were  found  to  be 
dependent  on  the  sizes  of  the  receiving  probes. 
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6"  cylinder  target 
3"  receiving  probe 


6"  cylinder  target 
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The  impulse  responses  of  a  6"  thin  cylinder  extracted  from  its 
scattered  fields  measured  by  three  different  receiving  probes  with 
lengths  of  2",  3"  and  4  1/2".  The  scattered  fields  were  regularized 
by  convolving  with  a  smooth  pulse  before  deconvolution.  The 
extracted  impulse  responses  of  the  target  were  found  to  be  quite 
independent  of  the  sizes  of  the  receiving  probes. 
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